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1.  Introduction 


This  report  discusses  several  constitutive  models  for  the  viscoelastic  response  of  nearly 
incompressible  solids.  All  of  the  models  presented  here  are  properly  invariant,  which  implies 
that  they  are  valid  for  arbitrarily  large  material  rotations.  They  should  prove  useful  for  moderate 
to  large  deformations,  depending  on  the  particular  model.  For  infinitesimal  deformations,  these 
finite  deformation  models  reduce  to  the  classical  theory  of  isotropic,  linear  viscoelasticity  with  a 
purely  elastic  pressure-volume  response.  Thermodynamic  effects  are  neglected  here,  although 
we  intend  to  incorporate  temperature  dependence  in  a  subsequent  report. 

The  term  “nearly  incompressible”  refers  to  materials  whose  response  in  shear  is  much  softer  than 
in  (purely  volumetric)  compression.  A  more  precise  definition,  at  least  for  elastic  materials,  is 
that  the  bulk  modulus  is  several  orders  of  magnitude  larger  than  the  shear  modulus.  For 
viscoelastic  materials,  the  shear  modulus  is  a  time-dependent  function  (the  stress  relaxation 
function),  but  a  similar  criterion  can  be  applied:  the  material  is  nearly  incompressible  provided 
the  bulk  modulus  is  several  orders  of  magnitude  larger  than  the  instantaneous  elastic  shear 
modulus  (see  section  3.2).  For  such  materials,  volumetric  and  shear  strains  of  the  same 
magnitude  will  generate  a  pressure  that  is  several  orders  of  magnitude  larger  than  the  shear  stress; 
conversely,  pressure  and  shear  stresses  of  the  same  magnitude  will  generate  a  volumetric  strain 
that  is  several  orders  of  smaller  larger  than  the  shear  strain.  Examples  of  nearly  incompressible 
materials  include  rubber,  biological  materials  with  high  water  content,  and  tissue  simulants  such 
as  ballistic  gelatin.  For  such  materials,  the  pressure  depends  primarily  on  the  current  volumetric 
strain;  in  particular,  the  viscoelastic  contribution  to  the  pressure  may  be  neglected  to  a  first 
approximation. 

Thus,  the  models  described  here  focus  on  the  viscoelastic  nature  of  the  shear  stress.  Any  such 
model  requires  that  this  stress  depend  not  only  on  the  current  value  of  the  strain  but  also  on  the 
past  values,  that  is,  on  the  history  of  the  strain  up  to  the  current  time.  One  way  of  doing  this  is  by 
means  of  a  hereditary  integral,  as  in  the  linear  theory  of  viscoelasticity  for  infinitesimal  strains;  cf. 
Wineman  and  Rajogopal  (/).  This  is  the  approach  taken  here.  However,  for  finite  deformations 
some  care  must  be  taken  to  ensure  a  properly  invariant  formulation.  This  is  most  easily  done  by 
working  with  the  material  (referential)  description  of  stress  and  strain.  Hence,  the  models 
considered  here  are  expressed  in  terms  of  the  second  Piola-Kirchhoff  stress  tensor  S  and  the 
Green-Lagrange  strain  tensor  E.  The  deviatoric  part  of  the  Cauchy  stress  tensor,  which  measures 
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shear  stress  in  the  current  configuration,  may  then  be  obtained  from  S  by  standard  relations.  The 
necessary  background  material  from  continuum  mechanics  is  summarized  in  section  2. 

In  section  3,  we  consider  a  linear  hereditary  integral  analogous  to  that  used  in  infinitesimal  strain 
theory.  The  general  theory  is  laid  out  in  sections  3. 1-3.2,  and  the  response  to  some  simple  strain 
histories  is  discussed  in  section  3.3.  In  section  3.4,  we  consider  the  Prony  series  approximation 
to  the  stress  relaxation  function.  This  approximation  has  two  main  advantages.  It  allows  for 
easier  calibration  of  the  model  to  experimental  data  (a  topic  not  addressed  here),  and  it  allows  one 
to  derive  incremental  relations  for  the  stress.  Several  forms  of  these  incremental  relations  are 
derived.  Some  of  these  are  analogous  to  relations  in  the  literature;  others  appear  to  be  new. 

Their  equivalence  for  piecewise-linear  strain  histories  and  for  arbitrary  strain  histories  is 
discussed.  These  incremental  relations  are  also  used  to  derive  explicit,  closed-form  solutions  for 
the  stress  for  several  simple  strain  histories. 

The  linear  model  in  section  3  cannot  be  expected  to  adequately  represent  the  material  response  for 
arbitrarily  large  deformations.  In  section  4,  we  consider  a  nonlinear  generalization,  which  retains 
some  of  the  useful  features  of  the  linear  model.  The  general  theory  is  outlined  in  section  4.1,  and 
a  simple  special  case  is  discussed  in  section  4.2.  In  section  4.3,  we  introduce  the  Prony  series 
approximation  to  the  general,  nonlinear  model  and  derive  incremental  relations  analogous  to 
those  for  the  linear  model.  Some  simple  special  cases  of  this  Prony  series  approximation  are 
discussed  in  section  4.4. 

The  incremental  relations  derived  here  (for  either  the  linear  or  nonlinear  models)  can  be  easily 
implemented  in  Lagrangean  finite  element  codes.  They  eliminate  the  costly  evaluation  of  a 
hereditary  integral  at  every  time  step.  The  explicit,  closed-form  solutions  for  the  stress  (for 
simple  strain  histories)  are  useful  for  verifying  the  model  implementation.  Intended  applications 
of  this  work  include  terminal  ballistics  simulations  involving  biological  materials  or  ballistic 
gelatin,  where  the  high  strain-rates  dictate  that  rate-dependence  in  the  material  response  be  taken 
into  account. 


2.  Stress  and  Deformation  Tensors 


We  follow  (for  the  most  part)  the  notation  and  terminology  in  Truesdell  and  Noll  (2),  Gurtin  (3), 
Bowen  (4),  and  Holzapfel  (5).  Let  F  denote  the  deformation  gradient,  and  let  J  denote  the 
Jacobian  of  the  deformation: 

J  =  det  F  .  (1) 
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Let  R.  U,  and  V  denote  the  rotation  tensor  and  the  right  and  left  stretch  tensors  in  the  polar 
decomposition  of  F: 

F  =  RU  =  VR .  (2) 

Let  C  denote  the  right  Cauchy-Green  deformation  tensor,  and  let  E  denote  the  Green-Lagrange 
strain  tensor: 

C  :=  FtF  =  U2  =  I  +  2E,  E:=^(C-I),  (3) 

where  I  denotes  the  identity  tensor.  Then 

C~  1  =  F~1F~t  =  U~2  =  I-2E  +  0(E2) ,  (4) 

where  F  T  denotes  the  inverse  transpose  of  F.  and 

J  =  V det  C  =  1  +  tr  E  +  0{E2) .  (5) 


Let  a  denote  the  Cauchy  stress  tensor,  <rn  the  rotated  (Cauchy)  stress  tensor,  P  the  1st 
Piola-Kirchhoff  stress  tensor,  and  S  the  2nd  Piola-Kirchhoff  stress  tensor.  They  are  related  as 
follows: 

crR  :=  RtctR  =  -U SU  =  -RTPU  ,  (6) 

u  U 

S  :=  JF-1crF-T  =  JU  ^ctrU  ^  =  F  lP  ,  (7) 

P  :=  JoF-J  =FS  =  JRctrU-1  ,  (8) 

<r  =  -FSFt  =  Rcr  Rt  =  -PFt  .  (9) 

J  R  J 

The  deviatoric  part  of  any  2nd  order  tensor  T  is  denoted  by  dev  T.  In  particular,  the  deviatoric 
stress  tensor,  dev  rr.  is  the  deviatoric  part  of  the  Cauchy  stress  tensor  and  represents  a  tensorial 
measure  of  shear  stress.  The  pressure  p  and  deviatoric  stress  tensor  are  given  in  terms  of  the 
Cauchy  stress  tensor  by 

p  \=  —  |trcr  and  dev  er  :=  cr  —  |  tr(er)J ,  (10) 

which  yields  the  standard  decomposition  of  the  Cauchy  stress  tensor  into  spherical  and  deviatoric 
parts: 

<r  =  —pi  +  dev  cr  .  (11) 
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On  substituting  this  into  equation  61,  we  see  that  the  rotated  stress  tensor  aR  can  be  expressed  in 
terms  of  the  pressure  and  deviatoric  stress  by 

crR  =  —  pi  +  devcrR  ,  deverR  =  RT(devcr)R.  (12) 

Conversely,  we  may  solve  equation  12  for  the  pressure  and  deviatoric  stress  in  terms  of  crR: 

p=—-trcrR,  dev  a  =  i?(dev  ctr)Rt  .  (13) 

o 

Similarly,  on  substituting  equation  11  into  equation  7i,  we  see  that  the  2nd  Piola-Kirchhoff  stress 
tensor  S  can  be  expressed  in  terms  of  the  pressure  and  deviatoric  stress  by 

S  =  -JpC~1  +  S*,  S*  :=  JF-\devcr)F-T .  (14) 

Conversely,  on  substituting  equation  9i  into  equation  10i  we  see  that  the  pressure  can  be 
expressed  in  terms  of  S  by 

P  =  --f  tr  ( FSFt )  =  ~  tr  (USU)  =  ~jj  tr  (SC)  .  (15) 

And  from  equation  14,  we  see  that  the  deviatoric  stress  can  be  expressed  in  terms  of  S*  by 

devcr  =  -FS*FT  =  - R(US*U)RT  .  (16) 

U  u 

It  follows  that  FS*Ft  and  US*U  are  deviatoric  tensors;  this  property  depends  on  the  fact  that 
the  scalar  p  in  the  decomposition  14i  of  S  is  the  pressure,  i.e.,  that  p  —  —  |  tr  cr. 

Finally,  we  note  that  there  are  relations  involving  the  1st  Piola-Kirchhoff  stress  tensor  P 
analogous  to  those  for  S  above,  but  we  do  not  list  them  here. 
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3.  A  Linear  Viscoelastic  Model 


3.1  General  Remarks 


We  consider  a  simple,  isotropic,  viscoelastic  constitutive  relation  expressed  in  terms  of  the  2nd 
Piola-Kirchhoff  stress  tensor  S  and  the  Green-Lagrange  strain  tensor  E;  thermodynamic  effects 
are  neglected  here.  The  advantage  of  working  with  these  stress  and  strain  tensors  is  that  any 
constitutive  relation  that  expresses  S(t)  as  a  function  of  the  history  of  E  (or  C  or  U )  up  to  time  t 
is  necessarily  properly  invariant;  cf.  (2). 


We  decompose  S  as 


S  =  -JpC~1  +  S, 


(17) 


and  assume  that  S  is  governed  by  the  isotropic,  linear  viscoelastic  constitutive  relation  25  below. 
The  decomposition  17  is  formally  similar  to  equation  lAi  but  is  not  equivalent  to  it.  Indeed,  the 
constitutive  relation  for  S  is  such  that  FSFT  and  U SU  will  generally  not  be  deviatoric  tensors, 
so  by  the  remark  following  equation  16  it  follows  that  the  scalar  p  in  equation  17  will  be  related 
to,  but  not  necessarily  equal  to,  the  pressure  p;  cf.  equations  23-24  below.  The  constitutive 
relation  for  p  will  be  chosen  in  such  a  way  that  the  pressure  p  depends  (either  exactly  or 
approximately)  on  the  current  volumetric  strain  only. 


The  other  stress  tensors  are  easily  computed  in  terms  of  p  and  S.  On  substituting  equation  17 
into  equation  82,  we  see  that 

P  =  -JpF-T  +  FS.  (18) 

Similarly,  on  substituting  equation  17  into  equations  9i  and  62,  we  see  that 


-pi  +  cr  , 

ir  :=  -  FSFT , 

u 

(19) 

-pi  +  aR  , 

aR  :=  jUSU. 

(20) 

where  neither  cr  nor  crR  are  deviatoric,  in  general. 
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The  above  relations  imply 


dev  cr  =  dev  cr  =  dev  ( FSFT )  , 
dev  crR  =  dev  crR  —  —  dev  (U  SU )  . 

u 


(21) 


(22) 


On  substituting  equations  19  and  20  into  equations  10i  and  13i,  respectively,  or  on  substituting 
equation  17  into  equation  15,  we  see  that  the  pressure  p  is  given  in  terms  of  p  and  S  by 


(23) 


where 


-  tr  (FSFt)  =  -  tr  (USU)  =  -  tr  (SC)  . 

U  u  U 


(24) 


tr  cr  =  tier R 


Equations  21-24,  together  with  the  spherical-deviatoric  decompositions  11  and  121?  give 
alternative  expressions  for  the  Cauchy  stress  tensor  cr  and  its  rotated  counterpart  crR. 

3.2  Constitutive  Assumptions 

3.2.1  The  Constitutive  Relation  for  S 

We  assume  that  the  value  of  S  at  time  t  is  given  in  terms  of  the  history  of  the  strain  tensor  E  up  to 
time  t  by  the  linear  single  integral  law 


S(t )  =  2G(Q)E(t)  + 


L 


— oo 


2 G(t-s)E(s)  ds 


(25) 


Here  G  :  [0,  oo)  — >  (0,  oo)  is  the  stress  relaxation  function  in  shear,  and  G  denotes  the  derivative 


of  G: 


GW  :=4G(s)> 

ds 


(26) 


so  that 


r\  r\ 

G(t  -s)  =  -G(t-s)  =  ~G( t  -  s ) . 


(27) 
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The  stress  relaxation  function  G  is  assumed  to  be  smooth,  convex,  and  strictly  decreasing  with  a 
positive  limit  at  infinity: 

G( 0)  >  G(s)  >  G( oo)  :=  lim  G(s )  >  0  ,  Vs  >  0  .  (28) 

s— >oo 

G( 0)  and  G( oo)  are  the  instantaneous  and  equilibrium  elastic  shear  shear  moduli ,  respectively. 

G  is  negative  and  strictly  increasing  and  asymptotes  to  0: 1 

G(0)  <  G(s)  <  G(oo)  =  0,  Vs  >  0  .  (29) 


For  the  remainder  of  section  3,  we  assume  that  the  material  is  undeformed  for  times  t  <  0,  so  that 


F{t)  =  R(t)  =  U(t )  =  C{t)  =  I  and  E(t)  =  0  for  t<  0. 


(30) 


Then  S(t)  is  zero  for  t  <  0,  and  the  constitutive  relation  25  for  S  reduces  to 


S(t)  =  2G(0)E(t)+  /  2 G(t-s)E(s)ds 


=  2G(0)E(t)  +  /  2 G(s)E(t-s)ds 

Jo 


(31) 


3.2.2  The  Constitutive  Relation  for  the  Pressure 

We  are  interested  in  materials  for  which  the  pressure  p  depends  primarily  on  the  current 
volumetric  strain.  To  impose  this  condition  exactly,  we  would  assume  a  (possibly  nonlinear) 
elastic  constitutive  relation  for  p  of  the  form 

P  =  P(J),  P(1)  =  0,  (32) 


where  the  condition  on  the  right  ensures  that  the  pressure  is  zero  in  the  undeformed  state,  where 


J  =  1.  Then 


P  =  «o 


(1- J)  +  C((l- J)2) 


(33) 


'Other  desirable  properties  of  G,  such  as  positive  relaxation  spectrum,  are  not  discussed  here  since  we  eventually 
assume  that  G  is  given  by  the  Prony  series  in  equation  64. 
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where  the  (initial)  bulk  modulus  k0  is  given  by 


K  0  =  “ 


dp 

dJ 


j= 1 


-P'(l)  • 


(34) 


Given  equation  32,  it  follows  from  equations  23-24  that  the  scalar  p  in  the  decomposition  17  of  S 
must  satisfy 

P  =  P( J)  +  ^ tr  ^  =  P(J)  +  tr  (SC)  .  (35) 

In  particular,  the  expression  on  the  right  introduces  a  viscoelastic  contribution  to  p.  With  S 
determined  from  equation  31  and  p  from  equation  35,  the  stress  tensors  P,  er,  and  crR  can  be 
determined  from  equations  18-20. 

If  one  is  ultimately  only  interested  in  computing  cr  or  crR,  then  there  would  be  no  need  to  actually 
compute  p:  simply  use  the  spherical-deviatoric  decomposition  11  or  12i,  with  p  determined  from 
equation  32,  and  dev  cr  or  dev  <jr  determined  from  S  via  equation  21  or  22. 2 

Next,  suppose  we  replace  the  constitutive  assumption  32  for  p  by  an  analogous  assumption  for  p: 

P  =  P(J ),  P(!)  =  0,  (36) 


so  that 


P  =  k  o 


(1  -  ./)  +  C((l  -  J)2 3) 


) 


with  k0  again  given  by  — p'(l).  Then  by  equations  23-24,  the  pressure  p  satisfies 


p  =  p{J) 


3^7  fr  (sc) 


(37) 


(38) 


It  is  clear  that  in  this  case  there  will  be  a  viscoelastic  contribution  to  the  pressure  (via  S ).  While 
such  an  effect  is  not  unreasonable  for  polymers,  there  is  no  reason  to  expect  that  it  would  be 
determined  by  the  stress  relaxation  function  G  in  shear  that  characterizes  S?  On  the  other  hand, 
for  nearly  incompressible  materials,  that  is,  materials  for  which  the  bulk  modulus  k0  is  several 
orders  of  magnitude  larger  than  the  instantaneous  elastic  shear  modulus  G(0),  the  trace  term  in 
equation  38  should  be  relatively  small  in  many  cases  of  interest,  so  that  p  m  p  ( ./) ,  i.e.,  we  recover 

2Of  course,  it  is  clear  from  equations  21  and  22  that  the  calculation  of  these  deviatoric  tensors  from  S  involves  the 
same  trace  calculation  as  in  equation  35,  so  this  approach  to  computing  cr  or  crR  does  not  involve  any  less  work  than 
the  one  in  the  preceding  paragraph. 

3One  reason  for  considering  equation  36  (in  conjunction  with  equations  17  and  31)  is  that  it  is  easy  to  construct 
a  viscoelastic  strain  energy  function  for  the  2nd  Piola-Kirchhoff  stress  tensor  S  in  this  case.  It  is  not  clear  that  this 
conclusion  necessarily  holds  when  p  is  given  by  equation  35,  i.e.,  when  equation  32  holds. 
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equations  32-34  as  an  approximation.  In  this  case  p  ~  p,  so  by  equations  14i  and  17,  we  see  that 
S  «  S*  also. 

The  constitutive  relation  17,  with  S  given  by  equation  31  and  p  given  by  either  equation  35  or 
equation  36,  is  isotropic  and  properly  invariant,  with  the  reference  configuration  a  natural  state  for 
the  material;  it  yields  the  elastic  pressure-volume  relation  p  =  p(  J),  either  exactly  or 
approximately.  For  terminal  ballistics  applications  with  nearly  incompressible  materials  like 
ballistic  gelatin,  for  which  the  bulk  modulus  is  four  orders  of  magnitude  larger  than  the  shear 
modulus,  this  modeling  framework  should  be  adequate  for  arbitrarily  large  rotations  and  for 
moderate  strains.  Whether  the  linear  constitutive  relation  31  for  S  suffices  to  characterize  the 
material  response  for  “large”  shear  strains  would,  of  course,  depend  on  the  particular  material. 
Some  nonlinear  constitutive  relations  are  discussed  in  section  4. 

Finally,  we  note  that  for  infinitesimal  deformations,  the  finite  deformation  theory  above  reduces  to 
the  classical  theory  of  isotropic  linear  viscoelasticity  with  a  purely  elastic  pressure-volume 
response. 

3.2.3  Jump  Discontinuities  in  Stress  and  Strain 

We  assume  that  the  strain  history  1 1— >  E(t)  is  piecewise  continuous.  Then  the  (equivalent) 
hereditary  integrals  in  equation  31  are  well-defined  and  are  continuous  functions  of  the  time  t. 

The  one-sided  limits  of  E  at  time  t  are  denoted  by 

E(t+)  \=  lim^(r)  and  E{t~)  :=  lim  E(t)  ,  (39) 

Tit  T^t 

and  the  jump  in  E  at  time  t  is  defined  by 

E^l(t)  :=  E(t+)  -  E(t~) ;  (40) 

similar  conventions  hold  for  the  other  variables.  Of  course,  if  E  is  continuous  at  the  instant  t 
then  E(t+ )  =  E(t~ )  =  E(t)  and  p£](t)  =  0.  But  if  E  suffers  a  jump  discontinuity  at  time  t, 
then  E(t)  need  not  be  defined. 

If  E  suffers  a  jump  discontinuity  at  time  t,  then  by  taking  the  one-sided  limits  of  equation  31  and 
using  the  continuity  of  the  integrals,  we  see  that  equation  31  holds  with  the  S(t)  and  E(t)  terms 
replaced  by  either  of  their  one-sided  limits.  On  taking  the  difference  of  these  limiting  relations, 
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we  see  that  the  integrals  cancel,  so  that  the  jump  in  S  at  time  t  is  given  by 


[S](f)  =  2G(0)  [£](*).  (41) 

In  particular,  past  jump  discontinuities  in  strain  do  not  affect  the  current  value  of  the  stress. 

The  condition  30  implies  that 

22(0")  =  0  and  5(0")  =  0 .  (42) 

We  assume  that  22(0)  equals  its  limiting  value  as  t  approaches  0  from  above.  Then  by 
equations  39-42  it  follows  that 

E(0)  =  E(0+)  =  122]  (0),  (43) 

5(0)  =  5(0+)  =  [51(0)  =  2G(0)22(0) .  (44) 


On  taking  the  rate  of  equation  31i,  one  finds  that  if  E  is  continuous  at  the  instant  t  but  its  rate 

•  _ 

(i.e.,  material  time  derivative)  E  suffers  a  jump  discontinuity,  then  S  is  also  continuous  at  time  t, 
but  its  rate  suffers  a  jump  discontinuity  given  by 

[S](t)  =  2G(0)  |[i](t) .  (45) 

In  view  of  equations  41  and  45,  we  see  that  jump  discontinuities  in  S  and  its  rate  are  governed  by 
the  instantaneous  elastic  shear  shear  modulus  G(0). 

3.2.4  Rate  Form  of  the  Relation  for  S 

To  simplify  things  further,  for  the  remainder  of  section  3  we  will  assume  that  the  deformation 
gradient  is  continuous,  except  possibly  for  a  jump  discontinuity  at  t  =  0;  then  the  same  holds  for 
the  strain  tensor  E.  This  assumption  does  not  lead  to  any  simplifications  of  the  viscoelastic 
constitutive  relation  31,  but  it  is  used  in  deriving  the  alternative  relation  47  below  in  terms  of  the 
strain-rate,  and  it  simplifies  the  incremental  forms  of  the  constitutive  relation  for  S  discussed  in 
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section  3.4.  Summarizing  the  assumptions  on  the  strain  history,  we  have 


Continuity  Assumptions  on  the  Strain  History: 


E(t)  =0  for  t  <  0  ; 


E(t)  is  a  continuous  function  of  t  for  t  >  0  ; 
E(t )  may  have  a  jump  discontinuity  at  t  —  0  ; 
E(t )  is  a  piecewise  continuous  function  of  t . 


(46) 


On  using  the  identity  27  in  equation  31  and  integrating  by  parts,  we  obtain  a  constitutive  relation 


for  S  in  terms  of  the  history  of  the  strain-rate  E: 


(47) 


For  continuous  strain  histories,  we  must  have  E(0)  =  0,  in  which  case  the  first  group  of  terms  on 
right  in  equation  47  drops  out: 


(48) 


This  should  be  appropriate  for  implementation  in  hydrocodes,  where  the  use  of  artificial  viscosity 
yields  continuous  strain  histories.  However,  for  the  discussion  of  certain  theoretical  results,  such 
as  stress  relaxation  tests  or  shock  waves  treated  as  propagating  jump  discontinuities,  we  need  to 
retain  the  possibility  of  an  initial  jump  in  strain.4 

4In  view  of  the  assumptions  46  on  the  strain  history,  which  are  essential  for  the  validity  of  equation  47,  it  follows 
that  this  equation  would  only  be  valid  for  a  shock  wave  arriving  at  time  0.  Note  that  no  such  restriction  applies  to  the 
original  relation  31.  If  strain  jumps  at  times  greater  than  zero  are  also  allowed,  then  the  relation  47  is  no  longer  valid: 
one  must  add  terms  involving  all  strain  jumps  up  to  time  f;  cf.  Scheidler  (6)  for  details. 
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3.3  Response  to  Special  Strain  Histories 
3.3.1  Stress  Relaxation  Tests 

Consider  an  ideal  stress  relaxation  test: 


F(t) 


I,  t  <  0  , 
Fq,  t  >  0 , 


(49) 


where  is  a  constant  tensor  (with  positive  determinant  Jo);  similarly  for  C,  R,  U,  and  V. 
Then 


0,  t  <  0, 

E(t)  =  {  (50) 

1  Eq,  t  >  0; 


so  by  equations  47  and  32  it  follows  that  in  a  stress  relaxation  test 


and 


__  [0,  t<  0, 

S{t)  = 

{  2G(t)E0,  t  >  0; 


lim  S(t)  =  2G(oo)  E0  ; 


pit) 


0,  t  <  0; 
p(J0),  t  >  0. 


(51) 

(52) 


(53) 


3.3.2  Constant  Strain-Rate  and  Strain-Rate  Jump  Tests 


Consider  a  continuous  strain  history  which,  for  some  time  f*  >  0,  satisfies 


m 


^1, 

0  <  t  <  t: 

e2, 

t  ^  t*  5 

(54) 


•  •  •  • 

where  E\  and  E-2  are  constant  symmetric  tensors.  If  E2  /  E\  then  there  is  a  strain-rate  jump  at 

time  f*.  The  strain  history  corresponding  to  equation  54  is  the  continuous,  piecewise-linear 
function: 

f  tE\,  0  <  t  <  U; 

m  =  {  .  .  (55) 

I  t*Ei  +  (t  —  t*)E 2,  t  > 
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When  Ei  and  E2  are  coaxial,  the  principal  stretches  (i.e.,  eigenvalues  of  U)  are  easily  computed 
from  equation  55,  but  they  are  not  piecewise-linear. 


By  equations  54  and  48,  we  have 


Sit )  = 


2 G(s)  ds 


Lao 


E  i, 


0  <  t  <  f* 


r  t-U 


2 G(s)  ds 


E  2  + 


2 G(s)  ds 


t  —  t* 


Ei,  f  >  f* 


(56) 


Hence  the  stress-rate  is  given  by 


2  G{t)Eu 


0  <  t  <  f*; 


S(t)  = 


2 G(f  —  QE2  +  2[G(f)  —  G(f  —  t*)]Ei,  t  >  f* 


Since  G(t)  asymptotes  to  G(oo )  as  t  — >  oo, 


(57) 


5(f)  — >  2G(oo)E2  ,  as  t  ^  oo 


(58) 


Note  that  for  times  f  such  that  0  <  f  <  t*  (i.e.,  the  top  lines  in  equations  54-57),  we  have  a 

•  • 

constant  strain-rate  test.  Of  course,  we  have  a  constant  strain-rate  test  for  all  times  if  E2  —  E i. 

For  the  special  case  where  E 2  =  0,  equation  55  reduces  to  a  linear  ramp  in  strain  over  the  time 
interval  [0,  f*]  to  a  fixed  strain  E *  :=  E( f*)  =  t*Ei. 

[  tEi,  0  <  f  <  f*; 

£(0  =  <  .  (59) 

I  t*Ei  =  E *,  f  >  f*. 


Since  =  0,  the  relation  562  for  5  (f )  for  times  £  >  £*  reduces  to 


5(f)  = 


2G(s)  ds 


t—t* 


El  =  r 


2G(s)  ds 


ft— £* 


i?* ,  f  >  f* . 


(60) 


For  real  stress  relaxation  tests  conducted  in  a  laboratory  (as  opposed  to  the  ideal  case  considered 
in  section  3.3.1),  the  strain  ramps  up  continuously  from  0  to  its  fixed  value  Et  over  some  short 
time  interval  [0,  f*].  The  strain  history  59  and  the  relation  60  for  5(f)  for  times  f  >  f*  provide 
better  approximations  to  such  tests.  Note  that  the  coefficient  of  Ei:  in  the  relation  on  the  right  in 
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equation  60  is  the  mean  value  of  G  over  the  time  interval  [t  —  f*,  f];  it  approaches  G(t)  as  A  — >  0, 
that  is,  as  the  ramp  time  approaches  zero.  Hence  for  fixed  f, 

5(f)  -►  2G(t)E*,  as  t*  — ►  0 ,  (61) 

which  is  consistent  with  the  relation  512  for  an  ideal  stress  relaxation  test  with  E{]  =  Ei:. 

3.4  Prony  Series  Approximation  and  Incremental  Relations 

In  section  3.4.1,  we  consider  the  Prony  series  approximation  to  the  stress  relaxation  function  G  ; 
cf.  equation  64.  When  this  is  substituted  into  the  hereditary  integrals,  5(f)  reduces  to  a  sum  of 
certain  simpler  hereditary  integrals.  This  sum  can  be  written  in  several  equivalent  ways,  three  of 
which  are  considered  here.  The  corresponding  simpler  hereditary  integrals  are  denoted  by 

An(t),  An(t),  and  An(f);  cf.  equations  68-75. 

The  strain  tensor  E(t)  can  be  computed  by  integrating  the  relation 

E  =  FtDF  ,  D  =  sym(V'u) ,  (62) 

where  Vv  is  the  spatial  velocity  gradient.  Of  course,  E(t)  can  also  be  computed  directly  from 
F(t)  using  equation  3,  with  F  obtained  by  integrating  the  relation  for  the  material  time  derivative 
of  the  deformation  gradient: 

F  =  (W)F.  (63) 

For  the  present  discussion,  the  strain  history  E  is  regarded  as  given  or  computed  at  the  times  of 
interest,  and  we  focus  attention  on  evaluating  the  integrals  *4"(f),  An(t),  and  A"(f).  Recall  that 
the  strain  history  satisfies  the  continuity  assumptions  46.  Sections  3.4.2  and  3.4.4  present  exact 
incremental  relations  for  An,  An,  and  A"  for  the  case  where  the  strain  history  is  piecewise-linear. 
In  section  3.4.5,  the  results  of  section  3.4.4  are  applied  to  some  special  piecewise-linear  strain 
histories.  In  section  3.4.6,  we  briefly  discuss  the  corresponding  approximate  incremental 
relations  for  arbitrary  strain  histories.  When  these  incremental  relations  are  substituted  into  the 
corresponding  expressions  for  5(f)  (cf.  equations  68,  71,  or  74),  we  obtain  incremental  relations 
for  the  stress.  Section  3.4.3  discusses  some  special  functions  that  appear  in  the  incremental 
relations. 
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3.4.1  Prony  Series  Approximation 

Henceforth,  we  assume  that  the  stress  relaxation  function  in  shear  is  given  (or  approximated)  by 
the  Prony  series 

N 

G(t )  =  /Too  +  fine  ;  (64) 

n=  1 

where  the  moduli  /iM  and  /in  and  the  relaxation  times  rn  are  positive.  Then 

N 

G(0)  /Too  T  ^  ^  fJ-n  •  A*0  )  G(oo)  /Too  ,  (65) 

n=l 


and 


Since  G(0) 


N 

G(t)  =  -J2~  e~t/Tn  ■ 

,  Tri 
n=l 

//o,  by  equation  44  we  see  that  the  value  of  S  at  time  zero  is  given  by 

S(0)  =  2n0E(0). 


(66) 


(67) 


On  substituting  equation  66  (with  t  — >  t  —  s)  into  the  constitutive  relation  31i  for  S(t)  and 
replacing  G(0)  with  //0,  we  obtain  the  following  relation  for  S(t )  at  any  time  t  >  0: 

N 

S(t )  =  2fj,0E(t)  -  J2  2^nAn(t) ,  (68) 

n=l 

where  the  dimensionless,  symmetric  tensor  A.n(t )  is  defined  by 

rt  i 

An(t):=  —e-^-s)/T-E{s)ds,  (69) 

JO 

so  that 

An(0)  =  0.  (70) 

Note  that  on  setting  t  =  0  in  equation  68  and  using  equation  70,  we  recover  equation  67. 

Similarly,  on  substituting  equation  64  into  the  constitutive  relation  47i  for  S(t),  we  obtain 

N  N 

S{t)  =  2/j,00E(t)  +  ^2 /inAn(t)  +  J22dne~t/TnE(0) ,  (71) 

n= 1  n= 1 
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where  the  dimensionless,  symmetric  tensor  An(t)  is  defined  by 


An(t) 


e-lt-s)h v.j5;(s)  ds, 


so  that 

An(0)  =  0. 

Note  that  on  setting  t  =  0  in  equation  71  and  and  using  equations  73  and  65,  we  recover 
equation  67. 

We  can  re-write  equation  71  as 


(72) 


(73) 


N 

S(t)  =  2fiooE(t)  +  2  AN  a n(t) ,  (74) 

n=  1 

where  the  dimensionless,  symmetric  tensor  A n(t)  is  defined  by 

An(A)  :=  An(t)  +  e~t/TnE(0) .  (75) 

Similarly,  on  substituting  the  expression  for  /i0  in  equation  65  into  equation  68  and  re-arranging, 
we  obtain  equation  74  with 

A n(t)  =  E{t) -An(t).  (76) 

It  is  also  easily  verified  directly  that  the  relations  76  and  75  for  A™  (A)  are  equivalent;  thus 

An(t)  +  e~t/TnE( 0)  =  E(t)  -  An(t) .  (77) 

This  may  be  used  to  express  An(t )  in  terms  of  An(t)  and  visa  versa;  then  equation  68  can  be 
obtained  from  equation  71  and  visa  versa. 

Given  that  G  is  represented  by  the  Prony  series  64  and  given  the  continuity  assumptions  46  on  the 
strain  history,  the  results  above  are  exact,  and  An  (t),  An(t),  and  A 71  (t)  are  continuous  functions 
of  t  for  t  >  0.  The  only  constitutive  property  that  they  depend  on  is  the  nth  relaxation  time  rn. 
Note  that 

An(0)  =  E{ 0) ;  (78) 

so  unlike  An  and  Arl,  the  tensor  An(0)  is  not  0  unless  E( 0)  =  0.  The  condition  E(0)  =  0  holds 
when  the  strain  history  is  continuous  for  all  times,  and  in  this  case  A n(t)  =  A" ( t )  by  equation  75, 
the  sum  on  the  far  right  in  equation  71  drops  out,  and  the  relations  71  and  74  are  indistinguishable. 
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In  the  constitutive  relations  68,  71,  and  74  for  S(t),  the  tensors  An,  An,  and  An  can  be  regarded 
as  internal  state  variables  that  modify  the  elastic  part  of  the  stress  (either  2 //0 77(7)  or  2/z0C_E(t)). 
One  can  easily  derive  evolution  equations  for  An  and  A"  by  differentiating  the  integral  relations 
69  and  72,  then  an  evolution  equation  for  An  can  be  obtained  from  equations  75  or  76.  However, 
for  purposes  of  computing  these  internal  state  variables,  it  is  simpler  to  work  directly  with  the 
integral  relations,  as  discussed  in  section  3.4.2. 

3.4.2  Incremental  Relations  for  Piecewise-Linear  Strain  Histories 

Here  and  in  section  3.4.4,  we  assume  that  in  addition  to  the  conditions  46,  E(t)  is  also  a 
piecewise-linear  function  of  time  t;  equivalently,  E  is  a  piecewise  constant  function.  Thus,  there 
is  some  sequence  of  times  tk  (k  =  0,1,2,...)  satisfying 

0  =  to  <  t\  <  t2  <  ■  ■  ■  <  tk- 1  <  tk  <  ■  ■  ■  (79) 

and  such  that  E  is  linear  (more  precisely,  affine)  over  each  time  interval  \tk- i,  tk].  We  set 

Aifc  =  tk  —  tk- 1 .  (80) 

Note  that  while  any  jump  discontinuities  in  the  strain-rate  E  must  occur  at  one  of  the  discrete 
times  tk,  E  need  not  suffer  a  jump  discontinuity  at  every  tk.  In  other  words,  E  may  be  linear 
(and  hence,  E  may  be  constant)  over  time  intervals  spanning  more  than  one  of  the  subintervals 
[tk- 1,  tk].  For  example,  for  the  strain-rate  jump  test  described  in  section  3.3.2,  we  need  only 
assume  that  tj  =  £*  for  some  j>  1  for  the  results  of  this  section  and  section  3.4.4  to  hold 
exactly.5  On  the  other  hand,  for  numerical  simulations  of  this  test  (or  of  more  general 
piecewise-linear  strain  histories)  with  the  tk  such  that  A tk  is  the  computational  time  step,  the  size 
of  A tk  would  be  limited  by  a  stability  condition,  so  that  many  subintervals  [tk- 1 .  tk\  might  be 
required  within  a  larger  time  interval  over  which  E  is  constant. 

We  claim  that  An,  An,  and  A™  are  given  exactly  at  the  times  in  equation  79  by  the  following 
incremental  relations  (here  k  =  1,2,...,  and  the  functions  a,  3,  and  /  below  are  defined  in 

5We  will  return  to  this  special  case  in  section  3.4.5. 
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section  3.4.3): 


An(tk)  =  e~At^An{tk- 1) 

+  a(Atk/rn)E(tk)  +  (3(Atk/rn)  E(tk_i) ,  (81) 

An(t0)  =  An(0)  =  0;  (81a) 

and 

An(tk)  =  e~Atk^An(tk_ i)  +  f(Atk/rn)  [. E(tk )  -  E(tk_ i)]  ,  (82) 

An(t0)  =  An(0)  =  0  ;  (82a) 

and 

A "(tk)  =  e-Atfc/r„An(4_i)  +  /(A4/rn)  [E(tk)  -  E( 4_i)]  ,  (83) 

A'1  (4)  =  An(0)  =  E(0) .  (83a) 


Note  that  the  incremental  relations  for  An  and  An  have  the  same  form  but  different  initial 
conditions. 

The  relations  81  and  82  are  derived  by  the  usual  procedure6  of  setting  t 
and  breaking  up  the  hereditary  integral  into  two  parts, 

rtk  rtk-i  rtk 

...  ds  —  ...ds+  . . .  ds  . 

J  0  Jo  J  tfc— i 

The  details  are  straightforward  and  not  included  here.7  The  relation  83  follows  from 
equations  82  and  75,  or  from  equations  81  and  76.  In  fact,  on  using  equations  75-77,  we  see  that 
any  one  of  the  incremental  relations  81-83  implies  the  other  two.  Now  set  t  —  tk  in  the  Prony 

6Cf.  Zocher  et  at.  (7)  and  Clements  (8).  Both  of  these  references  consider  the  more  general  case  where  the 
hereditary  integral  involves  a  reduced  time  variable. 

7But  the  following  observations  are  of  interest.  The  piecewise-linearity  assumption  is  not  needed  for  the  eval¬ 
uation  of  the  integral  over  [0,4-i];  this  integral  yields  first  term  on  the  right  in  equations  81  or  82  without  any 

assumptions  on  E  or  E.  However,  the  the  linearity  of  E  over  [4-i,ffc]  is  crucial  for  the  evaluation  of  the  second 
integral  on  the  right,  which  yields  the  bottom  right  group  of  terms  in  equation  81  and  the  second  group  of  terms  on  the 
right  in  equation  82. 


=  4  in  equations  69  or  72 

(84) 
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series  relations  68,  71,  and  74  for  S(t).  Then  for  k  —  0, 1,  2, ...  we  have 


N 


S(tk)  —  2/ro E{tk)  —  2/j,nAn(tk) , 

n— 1 
N 

5(4)  =  2/Xqo E{tk)  +  ^  ]  2 finAn  (4)  - 

n=  1 
N 

5(4)  =  2/i0OE(tk)  +  y^2/xwAw(4) . 


iV 


5^2/i, 


3  tk/^n 


72—1 


£(0), 


72=1 


(85) 

(86) 

(87) 


These  relations,  together  with  the  incremental  relations  81-83  for  A.n(tk ),  An(tk),  and  A" (4), 
yield  exact  (and  hence  equivalent)  relations  for  5(4)  for  piecewise-linear  strain  histories. 

Finally,  we  emphasize  that  that  the  continuity  assumptions  on  the  strain  history  (cf.  equation  46) 
have  been  used  in  the  deriving  the  incremental  relations  above  and  in  sections  3.4.4-3.4.6  below. 
Appropriate  modifications  when  these  conditions  are  relaxed  are  discussed  in  Scheidler  ( 6 ). 


3.4.3  The  functions  /,  a,  and  f3 


The  function  /  in  the  incremental  relations  82  and  83  for  An  and  An  is  defined  by 

1  -  e-* 

f{?)  = 


x 


x  ^  0  ; 


(88) 


1, 


x  —  0 . 


The  functions  a  and  (5  in  the  incremental  relation  81  for  An  are  defined  by 


a{x)  =  1  —  f{x)  and  (3(x)  =  f(x)  —  e  x  . 


(89) 


Then 


and 


a(x) 


e  x  +  x  —  1 
x 


x^O, 


0, 


x  =  0 ; 


P(x ) 


1  —  (x  +  l)e  x 
x 


x  ^  0, 


0, 


x  —  0 . 


(90) 


(91) 


See  figures  1  and  2.  /,  a,  and  f3  are  analytic  functions  on  the  real  line,  but  we  only  need  fix), 
a{x),  and  (3{x)  for  x  >  0,  and  these  all  lie  between  0  and  1.  /  is  strictly  decreasing  and 
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Figure  1.  f(x),  a(x),  and  (3{x)  for  0  <  x  <  20,  with  other  functions  shown  for 
comparison. 


asymptotes  to  0.  a  is  strictly  increasing  and  asymptotes  to  1.  3(x)  increases  from  0  at  x  — 
an  absolute  maximum  of  0.2984  . . .  at  x  —  1.793  . . .,  and  then  decreases  asymptotically  to  0. 
Indeed, 

P{x)~f(x) - ,  x  >  oo  , 

x 


0  to 


(92) 


and 


a(x)  ~  1 - ,  x 

x 


oo . 


(93) 


For  small  x,  we  have  the  approximations 


f{x)  =  1  -  ~x  +  \x2 - ,  (94a) 

2  o 

a(x)  —  \x  -  \x2  H - ,  (94b) 

2  o 

(3{x)  =  ]-x  -  ^ x 2  H - .  (94c) 

Zj  o 


3.4.4  Alternative  Relations  for  Piecewise-Linear  Strain  Histories 

In  this  and  the  next  paragraph,  we  assume  that  k  —  1,2,....  Given  the  assumptions  in 
section  3.4.2,  E  is  constant  over  any  time  interval  (tk-i,tk),  and  we  denoted  this  constant  value 
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f(x) 


0  1  2  3  4  5 


Figure  2.  f(x),  a(x),  and  /3(x)  for  0  <  x  <  5,  with  other  functions  shown  for 
comparison. 


•  •  • 

by  Ek.  However,  E  may  suffer  a  jump  discontinuity  at  any  of  the  times  tk,  so  12(4)  and 

•  • 

22(4- 1)  need  not  be  defined,  although  the  one-sided  limits  of  E(t)  as  t  —>  4  from  below  or  as 
t  — »  4-1  from  above  necessarily  exist  and  equal  Ek.  Clearly,  then,  we  have  the  following 
expressions  for  Ek\ 


E(tk)  '■=  lim-E(f), 

t  T  tk 

E(tU)  ■=  Ijm  E(t) , 

t i tfc-1 


E(t) ,  t  g  (4-ij  4)  — :  Ek  —  < 


E{t~k)  +  E(tU) 


E 


tk  +  tk- 1 


E(tk )  —  E(tk-i) 

tk  tk— 1 


E(tk )  —  E(tk- 1) 

Atk 


(95) 
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The  last  relation  above  yields 


E(tk)  —  E(tk- 1)  =  A tkEk  . 


(96) 


When  this  is  substituted  into  the  incremental  relation  82  for  An  or  the  incremental  relation  83  for 
A",  we  obtain  alternative  forms  of  these  relations.  In  particular,  equation  82  is  equivalent  to  the 
incremental  relation 


An{tk)  =  e~At^An{tk_ i)  +  r„(l  -  e~At^)Ek  ,  (97) 

An(t0)  =  An(0)  =  0,  (97a) 

where  Ek  is  given  by  any  of  the  expressions  in  equation  95. 8 

When  equation  97  is  substituted  into  equation  86,  we  obtain  another  (equivalent)  incremental 
relation  for  S(tk)  in  terms  of  E{tk)  and  An{tk).  This  relation  may  also  be  written  (for 
k  —  0, 1,  2, . . .)  as 

N 

S(tk )  =  2/j,00E(tk)  +  'y  ^  Sn(tk)  + 

n=l 

where  for  k  —  1,  2, . . .  we  have 

Sn(tk)  =  e~At^Sn(tk. i)  +  2/inrn(l  -  e~At^)Ek  ,  (99) 

and 


N 


y  ^  2// n 


t k/Tn 


n=  1 


^(0), 


(98) 


Sn(t0)  =  Sn(Q)  =  0.  (99a) 

Since  t0  =  0,  on  setting  k  =  0  in  equation  98  and  using  equations  99a  and  65,  we  recover 
equation  67.  When  the  strain  is  continuous  for  all  times,  so  that  22(0)  =  0,  the  incremental 
relation  98  reduces  (for  fc  =  0,l,2,...)to 


N 

S(tk)  =  2/j,00E(tk)  +  Y,  Sn(tk) ,  if  E(0)  =  0.  (100) 

77.— 1 


8Of  course,  equation  97  with  Ek  given  by  the  bottom  expression  in  equation  95  just  gives  back  the  relation  82. 
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3.4.5  Special  Piecewise-Linear  Strain  Histories 


Here  we  apply  the  results  at  the  end  of  the  preceding  section  to  some  special  piecewise-linear 
strain  histories. 


Constant  Strain-Rate:  We  assume  that  the  strain-rate  has  a  constant  value  E 1  for  t  >  0,  but 
allow  the  possibility  of  a  jump  in  strain  at  time  t  —  0,  so  that  -E'(O)  is  not  necessarily  0.  Thus 


E(t )  =  E i ,  t  >  0  , 


(101) 


and  the  strain  is  a  linear  function  of  time: 


E(t)  =  E(0)  +  tE1 ,  t>0. 


(102) 


Since  E  is  linear  over  the  time  interval  [0,  ti\  for  any  t  \  >  0,  and  since  At\  =  t\  —  t0  =  t\,  on 
setting  k  =  1  in  equation  99  and  using  equation  99a,  we  obtain 


Sn(ti )  =  2/inrn(l  -  e  tl/Tn)E1 . 


(103) 


Then  on  setting  k  —  1  in  equation  98  and  using  equations  103  and  102,  and  then  replacing  t\  with 
t,  we  see  that  for  any  t  >  0, 


S(t)  = 


N 


2^  +  E2^(l^  t/T" ) 


n= 1 


Ei  + 


N 


2 /Too  T  ^  ^  2  y  / , 


,-t/Tn 


n=L 


E{ 0) .  (104) 


This  also  holds  for  t  =  0,  since  it  reduces  to  equation  67. 


For  purposes  of  checking  the  implementation  of  the  incremental  algorithms  in  codes,  it  is  also 
useful  to  have  explicit  relations  for  A."  (t)  and  An(t)  in  this  special  case.  Proceeding  as  above 
and  using  equations  81-82  and  the  the  relations  for  /,  a,  and  (5  in  section  3.4.3,  we  find  that 


An(t)  =  a(t/rn)E(t )  +  P(t/rn)E(  0) 

=  rn(e-t/T-  +  t/rn  -  1  )E1  +  (1  -  e~t/Tn)E{ 0) , 


(105) 


and 


An(t)  =  rn(l  —  e~t^Tn)Ei 


(106) 


23 


Of  course,  these  relations  can  also  be  obtained  directly  from  the  definitions  69  and  72.  Also,  note 
that  on  substituting  equation  105  into  equation  68  and  using  equations  102  and  65,  we  recover  the 
relation  104  for  S(t).  Similarly,  we  can  recover  this  relation  by  substituting  equation  106  into 
equation  71  and  using  equation  102. 

For  the  case  where  there  is  no  initial  jump  in  strain  we  have  E(0)  =  0,  so  the  second  group  of 
terms  on  the  right  in  equations  104  and  105  drop  out.  The  solution  in  this  case  can  also  be 
obtained  by  substituting  the  Prony  series  64  into  the  corresponding  solution  56i  for  a  general 
stress  relaxation  function. 

• 

For  a  stress  relaxation  test  we  have  12(0)  =  E0  ^  0  and  E-v  =  0,  so  the  first  group  of  terms  on 
the  right  in  equations  104  and  1052  drop  out,  and  An(t)  =  0.  The  solution  in  this  case  can  also 
be  obtained  by  substituting  the  Prony  series  64  into  the  corresponding  solution  51  for  a  general 
stress  relaxation  function. 


Bilinear  Strain  (a  single  strain-rate  jump):  Here  we  assume  that  there  is  no  initial  jump  in 
strain,  so  that  12(0)  =  0,  and  that  the  strain-rate  satisfies 


E(t) 


E  i, 

0  <  t  <  t 

e2, 

c-F 

V 

c-F 

* 

(107) 


•  • 

so  that  there  is  a  strain-rate  jump  at  time  t*  if  E->  ^  E\.  This  is  the  case  considered  in 
section  3.3.2  (for  a  general  stress  relaxation  function  G).  Thus  the  piecewise-linear  strain  history 
is  given  by  equation  55. 


For  0  <  t  <  t * ,  the  strain-rate  is  constant  and  the  results  above  apply.  Thus  for  0  <  t  <  t * ,  S(t ) 
is  given  by  equation  104  with  12(0)  =  0.  It  remains  to  determine  S(t)  for  t  >  t*.  We  set 
1 1  =  t*.  Then  E  is  linear  over  the  time  interval  [t\,  i2]  for  any  t2  >  t*,  and  A t2  =  f2  —  £*.  On 
setting  k  =  2  in  equations  100  and  99,  and  then  replacing  t2  with  t,  we  see  that  for  any  t  >  t*, 


N 

S(t )  =  2/j,00E(t)  +  ^  [e-(*-**)A -»5n(i*)  +  2/inrn(l  -  e-(^*)^)i2]  , 

n=  1 


(108) 


where  Sn(t *)  is  given  by  equation  103  with  t\  =  t*.  On  substituting  that  relation  into 
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equation  108  and  using  the  relation  552  for  E(t),  we  find  that  for  any  t  >  t*, 


S{t)  = 


N 


2/Xqq f*  +  y^2/inrn(l  -  e  **/r")< 


,-(t-U)/Tn 


n= 1 


E  i 


JV 


2/Xoo(*  -  Q  +  ^2^nTn(l  “  e  **)/T") 


(109) 
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n= 1 

This  solution  can  also  be  obtained  by  substituting  the  Prony  series  64  into  the  corresponding 
solution  562  for  a  general  stress  relaxation  function. 


Similarly,  we  obtain  the  following  explicit  relations  for  A.n(t )  and  An(t )  for  t  >  i*: 


An(t)  =[t*  +  rn(e-'*/r"  - 


(HO) 


and 

An(t)  =  t„(1  -  e-u/Tn)e -(t~t*VT"E1  +  rn(l  -  e-^-u)/Tn)E2 .  (Ill) 

Note  that  on  substituting  equation  110  into  equation  68  and  using  equations  552  and  65,  we 
recover  the  relation  109  for  S(t).  Similarly,  we  recover  this  relation  by  substituting  equation  111 
into  equation  71  and  using  22(0)  =  0  and  equation  552. 

• 

For  the  special  case  where  E2  =  0,  we  have  a  linear  ramp  in  strain  over  the  time  interval  [0,  i*]  to 

•  _ 

a  fixed  strain  22*  :=  22 (i*)  =  t*E Then  E(t)  reduces  to  equation  59,  and  for  t  >  £*,  S(t)  is 

given  by  the  top  line  in  equation  109,  which  may  also  be  written  as 


S{t) 


N 


2/^oo  +  ^  ^  2  nnTn  (l 


e-U/Tn^e-{t-U)/Tn 


n= 1 


E * 


(112) 


These  formulas  can  also  be  obtained  by  substituting  the  Prony  series  64  into  the  corresponding 
solution  60  for  a  general  stress  relaxation  function. 

•  • 

Another  special  case  where  equation  109  simplifies  occurs  when  222  and  22,  differ  only  in  sign: 


22 2  —  —Ei . 


(113) 
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Then  for  t  >  £*,  equation  109  reduces  to 


S(t) 


N 

2/Xoo(2f*  ~t)  +  Y^  2 /^rn  (2e_(t_t*)/Tn  -  e_t/r’ 

n=l 


(114) 


3.4.6  Approximate  Incremental  Relations  for  General  Strain  Histories 

The  results  in  sections  3.4.2  and  3.4.4  are  exact  for  piecewise-linear  strain  histories,  and  there  is 
no  requirement  that  A tk  be  small.  Now,  we  drop  the  assumption  that  the  strain  history  is 
piecewise-linear  and  retain  only  the  continuity  assumptions  46.  The  above  results  then  give 
approximate  incremental  relations  for  A.11,  An,  and  An,  and  consequently  for  S.  provided  A tfc  is 
small  enough  that  E  is  approximately  linear  over  the  time  interval  [tk- 1,  tk]  ■  Here,  we  briefly 
discuss  these  approximations  (so  keep  in  mind  that  the  “=”  is  now  replaced  by  “ps”,  although  we 
won’t  bother  to  re-write  the  equations). 

The  incremental  relations  81-83  for  An,  An,  and  A”  are  equivalent  approximations.  That  is,  any 
one  of  these  approximations,  together  with  the  relations  75-77  and  equation  89,  implies  the  other 
two.  Consequently,  they  yield  equivalent  approximate  incremental  relations  for  S. 

• 

However,  since  E  is  not  necessarily  constant  over  the  time  interval  [tk-i,  tk],  the  expressions  on 
the  right-hand  side  of  equation  95  generally  yield  different  approximations  for  E(t)  on  this  time 
interval.  So  equation  97,  with  Ek  given  by  any  one  of  these  expressions,  generally  yields  a 
different  approximation  for  An(tk),  and  hence,  for  Sn(£fc).* 9  And  it  is  only  the  bottom  expression 
for  Ek  in  equation  95  which  yields  an  approximation  incremental  relation  that  is  equivalent  to 
those  in  section  3.4.2. 


4.  A  Nonlinear  Viscoelastic  Model 


We  continue  to  assume  a  decomposition  of  the  2nd  Piola-Kirchhoff  stress  tensor  S  of  the 
form  17,  with  constitutive  relations  for  p  and  the  pressure  p  as  discussed  in  section  3.2.2.  But 
now  we  consider  a  nonlinear  generalization  of  the  linear  viscoelastic  constitutive  relation  25  for 
S.  As  before,  the  Cauchy  stress  tensor  <r,  the  rotated  Cauchy  stress  tensor  crR,  and  their 
deviatoric  parts  can  then  be  obtained  from  equations  19-  22. 

•  •  • 

9The  relations  99-100,  with  the  top  expression  for  Ek  in  equation  95,  namely  Ek  =  E(tp),  are  analogous  to  the 
relations  in  Brad  Clements’  notes  (8)  for  the  special  case  where  the  reduced  time  there  is  the  ordinary  time. 
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The  nonlinear  constitutive  relation  for  S  is  introduced  in  section  4.1.1,  assuming  only  that  the 
strain  history  is  piecewise  continuous.  Stress  relaxation  tests  are  discussed  in  section  4.1.2.  An 
alternative  relation  for  the  stress  in  terms  of  the  history  of  E,  which  is  only  valid  when  the  strain 
history  satisfies  the  continuity  conditions  in  equation  46,  is  discussed  in  section  4.1.3.  The 
special  case  of  separable  strain  and  time  dependence  is  considered  in  section  4.2.1,  along  with  a 
brief  discussion  of  the  Prony  series  approximation  and  incremental  relations  for  this  special  case. 

In  section  4.3,  we  consider  the  Prony  series  approximation  to  the  general  nonlinear  model,  which 
offers  more  flexibility  than  the  special  case  considered  in  section  4.2.1,  and  discuss  the 
corresponding  approximate  incremental  relations.  Some  special  cases  of  this  general  Prony 
series  approximation  are  considered  in  section  4.4. 

Prior  to  section  4.3.2,  we  assume  only  that  the  strain  history  is  piecewise  continuous  unless 
specifically  stated  otherwise.  From  that  point  on,  we  re-introduce  the  continuity  and  smoothness 
conditions  46,  since  these  are  used  in  deriving  the  incremental  relations. 

4.1  The  General  Model 

4.1.1  Basic  Assumptions 


Consider  the  viscoelastic  constitutive  relation: 


S(t)  =  G(E(t),  0)  +  /  G(E(s),t- s)ds 

J  —OC 

=  G(E(t),  0)  +  /  G  (E(t  -  s),  s)  ds , 


(115) 


where 


G  :  Sym  x  [0,  oo)  — >  Sym 


(116) 


is  a  smooth  (and  generally  nonlinear)  function  of  both  arguments.  Here  Sym  denotes  the  space  of 
symmetric  tensors,  and  G  denotes  the  partial  derivative  of  G  with  respect  to  its  2nd  (temporal) 
argument: 

G(A,s)  :=^-G(A,s)  (117) 

for  any  (fixed)  symmetric  tensor  A.  We  assume  that  G  satisfies 


G(0,  s)  =  0  ,  Vs  >  0 ; 


(118) 
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then  also 


G(0,  s)  =  0  ,  Vs  >  0  . 


(119) 


The  conditions  118  and  119  guarantee  that  S(t)  =  0  for  t  <  t0  if  E(t )  =  0  for  t  <  t0,  consistent 
with  the  assumption  that  the  reference  configuration  is  a  natural  state.  Note  that  equations  115, 
17,  and  35  or  36  combine  to  give  a  properly  invariant  constitutive  relation.  A  viscoelastic 
constitutive  model  characterized  by  the  above  relations  is  often  referred  to  as  a  Pipkin-Rogers 
model 10  ;  cf.  (9).  The  material  is  isotropic  iff  G  is  an  isotropic  function  of  its  tensor  argument. 

The  linear,  isotropic  model  considered  in  section  3  is  a  special  case  of  the  above:  if 

G(E,s)  =  2G(s)E,  (120) 

then 

G{E,s)  =  2G(s)E,  (121) 

and  equation  115  reduces  to  the  constitutive  relation  25.  In  fact,  with  certain  qualifications, 

2 G(s) E  is  the  first-order  approximation  to  G (E.  s ).  We  have 

G(E,s)=  G(0,s)  +  DG{0,s)[E]  +  Os(E2),  (122) 

where  Os(E2)  denotes  a  function  of  E  and  s,  which,  for  fixed  s,  is  of  order  E2.  Here  DG{ 0,  s) 
denotes  the  derivative  of  G  with  respect  to  its  first  (tensor)  argument,  evaluated  at  0;  it  is  a 
fourth-order  tensor  that  acts  linearly  on  E  to  produce  a  second-order  tensor.  If  we  assume  that 
the  material  is  isotropic,  then 

DG{0,s)[E]  =  2G(s)E  +  X(s)(tiE)I 

for  some  functions  G  and  A.  But  since  tr  E  is  approximately  the  volumetric  strain  (cf. 
equation  5),  and  since  the  response  to  volumetric  strain  is  to  be  modeled  by  an  elastic  relation  for 
p  or  p  (cf.  section  3.2.2),  for  the  nearly  incompressible  materials  considered  here,  it  seems 
reasonable  to  restrict  attention  to  the  case  where  A  =  0.  On  using  this  assumption,  together  with 

ll)Typically,  this  term  is  applied  to  models  for  which  the  entire  2nd  Piola-Kirchhoff  stress  tensor  S  is  given  by  a 
relation  analogous  to  equation  115,  i.e.,  when  equation  115  holds  with  S  replaced  by  S  and  G  replaced  by  some 

j.  •  _ 

function  Q\  S(t)  =  Q  [E(t),  0)  +  f_oo  Q  ( E(s),t  —  s)  ds.  But  the  case  where  S  is  given  by  equation  17  with  S 

and  p  given  by  equations  115  and  36  can  be  obtained  from  this  by  taking  Q(E.  s')  =  —  Jp(J)C_1  +  G(E ,  s)  and 

•  • 

noting  that  Q(E,  s)  =  G(E,s). 
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the  condition  118,  we  see  that  equation  122  reduces  to 


G{E,s)  =  2G(s)E +  Os{E2) .  (123) 

Then  the  linear  constitutive  model  considered  in  section  3  can  be  regarded  as  a  first-order 
approximation  to  the  nonlinear  model  considered  in  this  section. 

Depending  on  the  particular  function  G,  it  may  be  simpler  to  express  G  (E,  s )  as  a  function  of  C 
and  s  using  the  relations  3,  say 

G(E,  s)  =  G(C,  s) .  (124) 

An  example  where  this  is  the  case  is  discussed  in  section  4.4.2.  Then 

G{E,s)  =  G{C,s) ,  (125) 

and  since  C  =  I  when  E  =  0,  the  conditions  118  and  119  imply  that 

G(I,s)  =  G(I,s)  =  0,  Vs  >  0 .  (126) 


Henceforth,  as  before  (cf.  equation  30),  we  assume  that  the  material  is  undeformed  for  times 
t  <  0,  so  that  E(t)  —  0  for  t  <  0.  Then  S(t)  =  0  for  t  <  0,  and  the  constitutive  relation  115 
reduces  to 

S(t)  =  G0(E(t ))  +  r 

Jo 

=  G0(E(t))  +  [‘ 

Jo 

where  we  have  introduced  the  instantaneous  elastic  response  function  Go,  defined  by 

Go(E)  :=  G{E ,  0)  =  G(C,  0)  =:  G0(C) ;  (128) 

in  view  of  equations  118  and  126  it  satisfies 

Go(0)  =  G0(J)  =  0.  (129) 


G ( E(s),t  —  s)  ds 
G  ( E(t  —  s) ,  s)  ds , 


(127) 


The  (equivalent)  hereditary  integrals  above  are  well-defined  for  strain  histories  that  are  only 
piecewise  continuous,  and  they  are  continuous  functions  of  t,  so  that  past  jump  discontinuities  in 
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strain  do  not  affect  the  current  stress.  If  E  has  a  jump  discontinuity  at  time  t,  then  equation  127 
holds  with  the  terms  S(t)  and  Eit.)  replaced  by  their  one-sided  limits.  Then  on  taking  the 
difference  of  these  limiting  relations,  we  see  that  the  integrals  cancel,  so  that  the  jump  in  S  at 
time  t  is  given  by 

(Sm  =  Go{E(t+))  -  G0(E(t~))  .  (130) 

As  before,  we  assume  that  E(0)  equals  its  limiting  value  as  t  approaches  0  from  above.  Then 
equation  43  holds,  and  equation  44  is  replaced  with 

5(0)  =  5(0+)  =  ffSj(O)  =  G0(E( 0))  ;  (131) 

that  is,  G0  (.E(O))  is  the  instantaneous  step  in  the  stress  S  corresponding  to  a  step  in  strain  of 
amount  E( 0)  from  the  undeformed  state. 

4.1.2  Stress  Relaxation  Tests 

Now  consider  the  stress  relaxation  test  50,  i.e.,  E(t)  =  E0  and  C(t)  =  C0  for  t  >  0  .  Then  for 
any  time  t  >  0, 

S ( t )  =  G (E0,  0 )  T  /  G (Eq,  t  —  s)  ds  ; 

Jo 

and  since 

G(Eo,  t  —  s)  =  —  —  G(Eq,  t  —  s) , 
ds 

the  integral  above  reduces  to  —  [G(_E0,  0)  —  G (E0,t)] .  Thus  for  any  time  t  >  0 ,  the  stress  S(t) 
in  a  stress  relaxation  test  is  given  by 

S(t)  =  G(E0,t)  =  G(C0,t).  (132) 

In  view  of  this  result,  G  and  G  are  called  the  stress  relaxation  functions  (cf.  (9)),  although  for  the 
linear  case  in  equation  120  it  is  more  common  to  refer  to  G  as  the  stress  relaxation  function. 
When  we  need  to  distinguish  between  these,  we  refer  to  G  as  the  scalar  stress  relaxation  function, 
and  to  G  and  G  as  the  tensor  stress  relaxation  functions.  Note  that,  in  principle,  the  tensor  stress 
relaxation  functions  are  completely  determined  by  the  response  in  ideal  stress  relaxation  tests.11 

By  setting  t  =  0  in  equation  132,  we  see  that  the  stress  at  time  zero  in  a  stress  relaxation  test  is 

"This  would  no  longer  be  the  case  for  a  more  general  single  integral  law  in  which  the  integrand  in  equation  127 
also  depends  on  the  current  value  E(t)  of  the  strain. 
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given  by  the  instantaneous  elastic  response  function: 


S{ 0)  =  G(E0, 0)  =  G0(E0)  =  Go(Co) ,  (133) 

consistent  with  equation  131. 

It  is  assumed  that  the  stress  S(t)  in  the  stress  relaxation  test  50  approaches  a  limiting  value 
S'(oo),  called  the  equilibrium  stress,  as  t  — >  oo.  Then  by  equation  132  it  follows  that 

S(oo)  =  lim  G(E0,t)  =  G(E0,  oo)  =:  Goo(^o) 

^°°  (134) 

=  G(C0,oo)  =:  Goo(G0). 

The  functions  G,x_  and  G^  defined  above  are  called  the  equilibrium  elastic  response  functions . 

By  equations  118  and  126,  we  have 


Goo(0)  =  Goo  (I)  =  0,  (135) 

which  implies  that  the  equilibrium  stress  corresponding  to  zero  strain  is  zero,  consistent  with  the 
assumption  that  the  reference  configuration  is  a  natural  state.  Note  that  the  equilibrium  and 
instantaneous  elastic  response  functions  are  related  by 

r°°  . 

Goo  (E0)  =  Gq(E0)  +  /  G(E0,s)ds ;  (136) 

Jo 

an  analogous  relation  holds  for  Goo,  G0,  and  G. 

4.1.3  Rate  Form  of  the  Relation  for  S 

The  constitutive  relations  in  section  4.1.1  are  valid  for  strain  histories  that  are  only  piecewise 
continuous.  For  any  instant  s  at  which  E(s)  exists,  the  chain  rule  gives 

^  G(E(s),t-s)  =  DG(E{s),t  -  s)[E(s)}  -  G(E(s),t-s )  .  (137) 

Recall  that  D  G  denotes  the  derivative  of  G  with  respect  to  its  first  (tensor)  argument,  so  that 
DG (E(s),t  —  s)  is  a  fourth-order  tensor  that  acts  linearly  on  the  strain  rate  E(s)  to  produce  a 
second-order  tensor. 

• 

On  solving  equation  137  for  G  (E(s),  t  —  s)  and  substituting  the  result  into  equation  127,  we  find 
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that 


S(t)  =  G(E(0),t)  +  [  DG(E{s),t-s)[E{s)}ds,  (138) 

Jo 

provided  the  strain  history  satisfies  the  stronger  conditions  in  equation  46.  In  particular,  for  the 
stress  relaxation  test  50,  for  which  22(0)  =  E 0  and  E(s)  =  0  for  s  >  0,  equation  138  reduces  to 
equation  132.  If  the  strain  history  is  continuous  for  all  times  (including  t  =  0),  then  we  must 
have  22(0)  =  0,  so  on  using  equation  118,  we  see  that  equation  138  reduces  to 

S(t)  =  [  DG(E(s),t- s)[E(s)]ds.  (139) 

Jo 

The  presence  of  DG,  rather  than  G,  in  the  integrands  above  may  result  in  a  more  complicated 
form  of  the  constitutive  relation  as  compared  with  equation  127.  An  exception  to  this  statement 
occurs  when  G  is  given  by  the  linear  relation  120;  then 

DG(E,s)[A]  =  2G(s)A  (140) 

for  any  symmetric  tensor  A,  and  equations  138  and  139  reduce  to  equations  47  and  48, 
respectively.  A  similar  simplification  in  equations  138  and  139  holds  for  the  separable  case 
considered  in  the  section  4.2. 

4.2  The  Separable  Case 

As  a  simple  example  of  the  general  theory  in  section  4.1,  we  consider  the  special  case  where  the 
stress  relaxation  function  G(E,  s )  is  a  separable  function  of  E  and  s. 

4.2.1  Constitutive  Assumptions  for  the  Separable  Case 

Assume  that  G  satisfies 

G(E,s)=2G{s)<S>{E),  (141) 

with  G  as  in  section  3,  and  $  :  Sym  — >  Sym  a  smooth  (and  generally  nonlinear),  dimensionless 
function  satisfying 

$(0)  =  0.  (142) 

The  material  is  isotropic  iff  is  an  isotropic  function.  In  this  case,  if  equation  123  also  holds, 
then  $  must  satisfy 

®(E)  =  E  +  0(E2) .  (143) 
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The  simple,  linear,  isotropic  model  considered  in  section  3  is  a  special  case  of  the  above  with  the 
remainder  term  in  equation  143  equal  to  0. 

From  equation  141  we  have 

G(E,s)  =  2G(s)&(E).  (144) 

And  from  equations  141,  128  and  134, 

G0(E)  =  2G(0)&(E),  G0O(E)  =  2G(oo)&(E).  (145) 


Then  the  general  constitutive  relation  127  reduces  to 

S(t)  =  2G(0)&(E(t))+  [  2 G(t-s)&(E{s))ds 

Jo 

rt 

=  2G(0)$(E(t))+  /  2G(s)&(E(t  —  s))  ds . 
Jo 

And  integration  by  parts  yields  the  alternative  form 

—  rl  d 

S(t)  =  2G(t)&(E(0))+  /  2 G(t-s)  —  &(E(  s))  ds , 

Jo  ds 


(146) 


(147) 


provided  the  strain  history  satisfies  the  stronger  conditions  in  equation  46.  When  22(0)  =  0,  this 
reduces  to 

S(t)=  /  2G(t-s)  —  &(E(s))  ds.  (148) 

Jo  ds 

The  relation  148,  or  something  analogous  to  it,  is  often  encountered  in  the  literature  on  nonlinear 
viscoelasticity. 


Note  that  equation  139  reduces  to  equation  148  when  G  is  given  by  equation  141.  Indeed,  in  this 
case,  we  have 

DG(E,s)[A]  =  2 G(s)D<f>(E)[A] ,  (149) 

for  any  symmetric  tensor  A,  so 


DC(£(s),f-s)[£(s)]  =  2G(t^s)D®(E(s))[E(s)] 


(150) 


4.2.2  Alternative  Notation  for  the  Separable  Case 


Depending  on  the  particular  form  of  the  function  <I>,  it  may  be  simpler  to  express  <I>(E)  as  a 
function  of  C  using  the  relations  3,  say  &(E)  =  Then  by  the  condition  142,  we  must 

have 

*(!)  =  0.  (151) 

Also,  for  use  in  the  relations  above,  as  well  as  for  the  discussion  below,  it  is  convenient  to 
introduce  the  notation 

£:=<&  (C)  =  ${E)  =  E  +  0{E2),  (152) 

where  the  approximation  on  the  right  is  just  equation  143.  Then  we  may  regard  £  as  a 
generalized  strain  tensor.  Since 


£(»)  =  *(«(»))  =  *(C(<0), 


(153) 


we  have 

£(»)  =  4S(B(s))  =  £.*(J5(S))[£(3)]  , 
and  similarly  in  terms  of  \Er.  With  this  notation,  equation  146  can  be  written  as 


S(t )  =  2G(0)£(t)  + 


2 G(t  —  s)£(s )  ds 


=  2G(Q)£{t)  + 


2 G(s)£(t  —  s)  ds; 


and  equation  147  can  be  written  as 


S{t )  =  2G(t)£{0)  + 
=  2G(t)£(0)  + 


2  G(t  —  s)£(s)  ds 
2  G(s)£(t  —  s)  ds; 


(154) 


(155) 


(156) 


and  then  equation  148  is  just  equation  156  with  5(0)  =  0.  These  relations  are  completely 
analogous  to  those  for  the  linear  theory  in  section  3;  cf.  equations  31,  47,  and  48,  respectively.  In 
particular,  when  <E »(E)  =  E,  so  that  £  =  E,  we  recover  the  relations  in  section  3. 
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4.2.3  Prony  Series  and  Incremental  Relations  for  the  Separable  Case 


Now  assume  the  Prony  series  approximation  64  for  G  in  equation  141.  Then 


G(E,t) 


N 


2  i-ioo  T  ^  ^  2  fir 


n= 1 


N 

=  2 fioo^E)  +  J22 line~t/Tn<f>{E) . 

n=  1 


(157) 


Also,  assume  the  continuity  conditions  46,  so  that  both  equations  155  and  156  hold.  Then  all  of 
the  relations  in  section  3.4.1  hold  with  E  replaced  by  £  everywhere.  Because  of  this  close 
correspondence,  we  will  not  bother  to  write  these  new  relations  here,  but  simply  refer  to  them  as 
the  “modified”  relations  in  the  discussion  below. 


Regarding  the  incremental  forms  of  the  modified  relations  for  AJ\  An,  and  A"  ,  first  note  that 
even  if  E(t)  is  a  piecewise  linear  function  of  t,  £{t )  generally  will  not  be,  due  to  the  possible 
nonlinearity  of  <&.  However,  if  E  is  approximately  linear  over  the  time  interval  [tfc-i,  A],  then  so 
is  £,  provided  At),  is  sufficiently  small.  When  this  is  the  case,  the  incremental  relations  81-83 
with  E  — >  £  yield  approximate  incremental  relations  for  the  modified  tensors  A" ,  An,  and  An. 
These  incremental  relations  are  equivalent:  any  one  of  them,  together  with  equations  75-77 
(with  E  — >  £)  and  equation  89,  implies  the  other  two.  Consequently,  they  yield  equivalent 
approximate  incremental  relations  for  S(t)  via  the  modified  form  of  the  relations  85-87  for 

S(tk).  Other  approximate  incremental  relations  can  be  obtained  from  those  in  section  3.4.4  via 

•  • 

the  replacements  E  — ►  £  and  E  —>  £. 

4.3  Prony  Series  Approximation  for  the  General  Model 

For  the  special  case  considered  in  section  4.2.3,  where  G  has  the  separable  form  141  with  G  given 
by  the  Prony  series  64,  the  stress  relaxation  function  is  given  exactly  by  equation  157.  This  can 
also  be  regarded  as  an  approximation  to  a  general  stress  relaxation  function  G  in  the 
Pipkin-Rogers  model.  Clearly,  a  better  approximation  could  be  obtained  if  we  allowed  the 
function  $  in  the  nth  term  of  equation  157  to  vary  with  n.  This  is  the  case  considered  here. 
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4.3.1  Prony  Series  Approximation  for  the  Stress  Relaxation  Function 


Assume  that  the  stress  relaxation  function  G  in  section  4.1.1  is  given  (or  approximated)  by  the 
Prony  series 

N 

G(E,t )  =  2/u00$00(.E)  +  fine~t/Tn^n{E) ,  (158) 

n— 1 

where  each  :  Sym  — >  Sym  (n  =  1, . . . ,  N,  oo)  is  a  smooth  (and  generally  nonlinear), 

dimensionless  function,  and  the  moduli  /i^  and  //„  and  the  relaxation  times  rn  are  positive.  Then 

N  9 

G{E,  t)  =  -J2~  e~t/Tn®n(E) .  (159) 

„_i  Tn 


The  material  is  isotropic  if  each  T»„  is  an  isotropic  function,  which  will  be  assumed  here. 
Depending  on  the  particular  form  of  the  function  <f>„,  it  may  be  simpler  to  express  &n(E)  as  a 
function  \Pn  of  C  using  the  relations  3;  an  example  where  this  is  the  case  is  discussed  in 
section  4.4.2.  Thus,  for  n  =  1, . . . ,  N,  oo,  we  have12 

*n(C)  :=  $n(E)  =  E  +  On{E 2) ,  (160) 


where  the  approximation  assumed  on  the  right  guarantees  that  the  approximation  123  holds  with 
G(s)  given  by  the  Prony  series  64: 


N 

G(E,t)  =  2G(t)E  +  Ot(E2),  G(t)=fj,0O  +  '%2(jLne-t/Tn.  (161) 

n=l 


In  particular, 


$n(0)  =  *„(!)  =0. 


(162) 


Thus,  as  the  strain  E  becomes  smaller,  the  Prony  series  approximation  158  for  the  general  stress 
relaxation  function  G  in  section  4.1.1  approaches  the  Prony  series  approximation  for  the  linear 
theory  considered  in  section  3.  It  follows  that  the  moduli  ^  and  jin  and  the  relaxation  times  rn 
can  be  determined  from  the  viscoelastic  response  to  small  deformations.  Departure  from  the 
linear  viscoelastic  theory  enters  only  through  the  time-independent  functions  T»„.  This  is  an 
advantage  (as  well  as  a  limitation)  of  the  approximation  158. 

uOn(E2)  denotes  a  function  of  E  and  n  which,  for  fixed  n,  is  of  order  E2. 
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From  equation  158  it  follows  that  the  equilibrium  elastic  response  function  is  given  by 


Goo(E)  :=  G(E,oo)  =  2ji00<&00(E) , 


(163) 


with  /loo  again  representing  the  equilibrium  elastic  shear  modulus.  We  call  the  normalized 
equilibrium  elastic  response  function.  The  instantaneous  elastic  response  function  is  given  by 

N 

Go(E)  G(E ,  0)  =  2/i00d>00(_E)  +  2p,nQ>n(E) .  (164) 

n— 1 


This  may  also  be  written  as 

Gq(E)  =  2/i0$0 (E) , 

where  the  instantaneous  elastic  shear  modulus  /i0  is  defined  as  in  equation  65,  i.e., 

N 

do  ■  doc  T  ^  ^  bn  j 

71=1 


(165) 


(166) 


and  the  dimensionless  function  d>0  is  defined  by 

N 

$0 (E)  :=  ^  ^oc(^)  +  J2  —  *n(E) .  (167) 

do  ^  do 

We  call  $o  the  normalized  instantaneous  elastic  response  function.  Note  that  the  sum  of  the 
dimensionless  coefficients  in  equation  167  is  unity,  so  on  using  the  small  strain  approximations 
from  equation  160  we  see  that  <E>0(_E)  =  E  +  O0(E2),  i.e.,  equation  160  also  holds  for  n  —  0. 
Finally,  note  that  in  view  of  equation  160i,  the  relations  158,  159,  163-165  and  167  hold  with  the 
replacements  G  —>  G,  E  —>  C  and  d?n  — ►  n  (n  =  0, 1, . . . ,  N,  oo). 

4.3.2  Relations  for  the  Stress 

On  substituting  equation  159  (with  E  — >  E(s)  and  t  — >  t  —  s)  and  equation  165  into  the 
constitutive  relation  127i  for  S(t),  we  obtain 


N 

S{t )  =  2/i0$o {E(t))  -J22 dnAn(t) , 

n=  1 

where  the  dimensionless,  symmetric  tensor  A.n(t )  is  defined  by 

rt  1 

An(t)  :=  /  —  e-{t-s)/Tn®n(E{s))  ds  , 

Jo  Tn 


(168) 


(169) 
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so  that 


An(0)  =  0 . 


(170) 


The  relations  168-170  generalize  the  relations  68-70  for  the  linear  case. 

Henceforth,  we  again  assume  the  continuity  and  smoothness  conditions  46  on  the  strain  history. 
Then  on  integrating  equation  169  by  parts,  we  obtain 

An(t)  =  &n(E{t))  -  e~t/T^n(E( 0))  -  An(t) ,  (171) 

where  the  dimensionless,  symmetric  tensor  An(t )  is  defined  by 

An(t )  :=  I  e-(‘-s)/r"^Un(£(S))  ds  ,  (172) 

so  that 

An(0)  =  0.  (173) 

On  substituting  equation  171  into  equation  168  and  using  equation  167,  we  obtain 

N  N 

S(t )  =  2/^4^))  +  ^2 HnAn{t)  +  J2‘2Vne-t/Tn®n{E(0))  .  (174) 

n=  1  72—1 

The  relations  174,  172,  and  173  generalize  the  relations  71-73  for  the  linear  case. 

4.3.3  Alternative  Notation 

The  expressions  above  for  the  tensors  S(t),  An  (/),  and  An(t),  as  well  as  the  approximate 
incremental  relations  for  them,  which  are  given  in  the  next  section,  can  be  simplified  somewhat 
by  using  introducing  the  following  notation:13 

Sn{s)  :=<i>n(E{s))=*n(C{s)),  n  =  0, 1, . . . ,  IV,  oo  .  (175) 

On  using  equation  175  in  equations  168  and  169,  we  obtain 

N 

S(t )  =  2fi0£\t)  -  2 »nAn(t) ,  (176) 

72—1 

/■*  i 

An(t)=  - e-{-t~s)lT-En(s)ds ,  (177) 

J0  Tn 

13Note  that  £n  does  not  denote  £  raised  to  the  nth  power. 
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which  more  closely  resemble  the  relations  68-69  for  the  linear  case. 

From  equation  175,  the  material  time  derivative  of  £ 71  is 

c*W==^"W  =  |*.(*W) 

=  D$n(£;(5))[i(s)]J 

and  similarly  in  terms  of  \Fn.  On  using  equations  178  and  175  in  equations  174  and  172,  we 
obtain 

N  N 

S(t)  =  2/j,00£0°(t)  +  J]2/rnAn(f)  +  ^  2^ne~t/T-£n{Q) , 

71=1  71=1 

rt 

An(t)  =  /  e-{t-s)/Tn£n(s)ds, 

Jo 

which  more  closely  resemble  the  relations  71-72  for  the  linear  case. 

Also,  by  equation  175  the  relation  171  between  A.n(t)  and  A"  (t)  can  be  written  as 

£n(t )  -  An(t)  =  An(t)  +  e~t/Tn£n(0)  =:  A n(t) ,  (181) 

which  more  closely  resembles  the  relations  75-77  for  the  linear  case.  On  using  the  relation  on 
the  right  above,  we  see  that  equation  179  can  be  written  as 

N 

S( t)  =  2n00£00(t)  +  2/inAn(f) ,  (182) 

n=  1 

analogous  to  the  relation  74  for  the  linear  case.  And  by  equations  181,  173  or  170,  and  175,  we 
see  that 

Ara(0)  =  £n(0)  =  &n(E{ 0))  =  ^n(C(0))  .  (183) 

The  condition  £)(0)  =  0  holds  when  the  strain  history  is  continuous  for  all  times,  and  in  this  case 
<&n(.E(0))  =  0  by  equation  162,  so  A"(0)  =  £n(0)  =  0  by  equation  183;  then  A n(t)  =  An(t )  by 
equation  181,  the  sum  on  the  far  right  in  equation  179  drops  out,  and  the  relations  179  and  182  are 
indistinguishable. 


(179) 

(180) 
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Finally,  note  that  the  definition  175  and  the  condition  160  imply 


en(t)  =  E{t)  +  On(E(t)2)  ,  n  —  0,1, . . . ,  N,  oo  . 


(184) 


Thus  we  may  think  of  the  £n  as  generalized  strain  tensors.  And  substitution  of  equation  184  into 
equations  176-177,  179-180,  and  181-182  yields  the  relations  for  the  linear  theory,  namely 
equations  68-69,  71-72,  and  74-77,  respectively,  as  small  strain  approximations. 

4.3.4  Approximate  Incremental  Relations 

Even  if  E(t)  is  a  piecewise  linear  function  of  t,  the  tensors  £n(t)  generally  will  not  be,  due  to  the 
possible  nonlinearity  of  the  functions  <F„  in  equation  175.  On  the  other  hand,  given  the 
continuity  conditions  46,  we  can  assume  that  the  sequence  of  times  4  in  equation  79  are  chosen 
so  that  not  only  E  but  each  £  n  is  approximately  linear  over  each  time  interval  [tk- 1,  4] .  Then, 
on  using  the  relations  for  A!1,  An,  and  A”  in  section  4.3.3,  and  proceeding  as  in  section  3.4.2,  we 
can  derive  approximate  incremental  relations  for  A1,  An,  and  A”.  By  equation  175,  we  have 


£n(tk)  =  &n(E{tk))  =  tfn(C(4)) 


(185) 


for  n  —  0, 1, . . . ,  N,  oo.  And  as  before,  we  let  A 4  =  tk  —  4-i-  Then,  we  find  that 


A\tk)  «  e-AWr„An(4_i) 

+  a(Atk/rn)£n(tk)  +  P(Atk/rn)£n(tk-i) , 


(186) 


An(t0)  =  A"(0)  =  0; 


(186a) 


and 


An(tk)  «  e-At^An(tk_ i)  +  f(Atk/rn)  [£n(tk)  -  £n(4-i)]  , 


(187) 


An(t0)  =  An(  0)  =  0; 


(187a) 


and 


An(4)  ~  e-Atfc/TnAn(4_i)  +  f(Atk/rn)  [ £n(tk )  -  £"(4-i)]  , 


(188) 


A"(4)  =  An(0)  =  £n(0) . 


(188a) 
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These  are  analogous  to  the  incremental  relations  81-83a  for  the  linear  case. 


Note  that  the  incremental  relations  for  An  and  A"  have  the  same  form  but  different  initial 
conditions.  Furthermore,  all  three  incremental  relations  are  equivalent:  any  one  of  them, 
together  with  the  relations  181  and  89,  implies  the  other  two.  Consequently,  they  yield 
equivalent  approximate  incremental  relations  for  S  via  equations  176,  179  and  182. 


Other  approximate  relations  may  be  obtained  as  in  section  3.4.4.  Since,  by  assumption,  £ n  is 

approximately  linear  over  each  time  interval  [tk-i,  tk],  its  rate  £n  will  be  approximately  constant 

•  • 

over  that  time  interval.  Let  (£n)k  denote  any  approximation  for  £"  over  [tk~ i,  tk\.  Five 
(generally  different)  approximations,  analogous  to  those  in  equation  95,  are  given  below: 


£n(4)  :  =  lim5n(f), 

t  I 

<£"*(<*-,)  :=  Mm  £”(*), 

4 i tfc-l 


tk  T  tk— 1 


£n(tk)  -  g"(4_i) 
tk  tk— 1 


£'\tk)  -  £n(tk- 1) 


(189) 


As  before,  we  have  allowed  for  the  possibility  that  £n  may  suffer  a  jump  discontinuity  at  tk 
and/or  tk-i.  The  last  choice  for  (£")k  above  is  equivalent  to  the  relation 

£n(tk)  -  £n(tk- 1)  =  A tk(£n)k  .  (190) 


When  this  is  substituted  into  the  incremental  relation  187  for  An  or  the  incremental  relation  188 
for  An,  we  obtain  alternative  forms  of  these  relations.  For  example,  equation  187  is  equivalent  to 
the  incremental  relation 

An(tk )  «  e~At^An(tk- 1)  +  rn(  1  -  e~At^)(£n)k  (191) 

when  ( £n)k  is  given  by  the  bottom  expression  in  equation  95;  this  is  analogous  to  equation  97. 
Then  equation  191  with  any  of  the  other  choices  for  (£n)k  in  equation  95  yields  a  (generally) 
different  approximate  incremental  relation  for  An.  These  relations  and  their  analogs  for  An  yield 
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approximate  incremental  relations  for  S  via  equations  179  and  182.  Finally,  one  can  proceed  as 
at  the  end  of  section  3.4.4  and  obtain  approximate  incremental  relations  for  S'  in  a  form  analogous 
to  equations  98-100. 

4.4  Special  Cases  of  the  Prony  Series  Approximation 

Regarding  the  Prony  series  approximation  158  for  the  stress  relaxation  function  G,  we  have  not 
placed  any  restrictions  on  the  functions  (n  —  1, . . . ,  N,  oo)  other  than  those  discussed  in  the 
first  paragraph  in  section  4.3.1.  In  particular,  the  relations  in  section  4.3  do  not  require  that  these 
functions  be  related  to  each  other  in  any  simple  way.  Here  we  consider  some  special  cases  of  this 
general  theory. 

4.4.1  A  Simple  Non-Separable  Case 

When  =  •••  =  <& at  =  <&oo  =:  <&,  the  Prony  series  158  reduces  to  equation  157,  and  we 
recover  the  Prony  series  approximation  to  the  separable  stress  relaxation  function  considered  in 
section  4.2.  A  special  case  of  equation  158  that  is  not  necessarily  separable  occurs  when 

<E»i  =  =  <E»at  =:  <f>  ,  (192) 


with  $  possibly  distinct  from  the  normalized  equilibrium  elastic  response  function  $00.  Then 
the  stress  relaxation  function  in  equation  158  reduces  to 


G(E,t)  —  2/x004?00(_E)  + 


&{E). 


As  before,  the  equilibrium  elastic  response  function  is  given  by  equation  163: 


(193) 


Goo(E)  :=  G(E,  00)  =  2/u00$00(£)  . 


(194) 


And  on  setting  t  =  0  in  equation  193  and  using  equation  166,  or  by  using  equations  192  and 
164-167,  we  see  that  the  instantaneous  elastic  response  function  is  given  by 

Go (E)  =  2/i00<&00(_E)  +  2(/i0  —  /ioo)*(E)  =  2/io4?o(-E0  ,  (195) 

where  the  normalized  instantaneous  elastic  response  function  is  given  by 

$0 (E)  =  <M£)  *(E) .  (196) 

^0  V  Fo  / 
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In  view  of  equation  192,  the  tensors  £n  defined  in  equation  175  are  identical  for  n  —  1, . . . ,  N: 


£n{t)  =  £{t)  :=  $ (E(t))  =  $ (i C(t ))  ,  n  =  l,...,N.  (197) 

However,  the  tensors  £°(t)  and  £°°(t),  which  appear  in  the  relations  176,  179,  and  182  for  S(t), 
generally  do  not  coincide  with  £(t ):  by  equations  175  and  196,  we  have 

£°°{t)  ■■=  &oo(E(t))  =  *oo(C(t))  ,  (198) 

and 

£\t)  :=  *0(E(t))=V0(C(t)) 

/  X  (199) 

Ho  \  Ho  ) 

When  $00  =  $  we  recover  the  separable  case  in  section  4.2.  And  when  $00(£;)  =  <&(E)  =  E, 
this  reduces  to  the  linear  case  in  section  3.. 


The  Simplest  Non-Separable  Case:  The  simplest  special  case  of  the  above,  which  does  not 
reduce  to  the  separable  case,  occurs  when  the  normalized  equilibrium  elastic  response  function 
$oo  is  allowed  to  be  nonlinear  but 

<$>(E)  =  E,  (200) 

in  which  case 

£n(t)  =  £(t)  =  E(t),  n  =  l  (201) 

In  other  words,  the  equilibrium  elastic  response  may  be  nonlinear,  but  the  viscoelastic  part  of  the 
response  is  linear  in  E. 

In  this  case  the  relations  177,  180,  and  181  for  A.n,  An,  and  An  reduce  to  the  relations  69,  72,  and 
75-76  of  the  linear  theory.  Hence  the  incremental  relations  for  A!\  An,  and  An  in  section  4.3.4 
reduce  to  those  of  the  linear  theory  as  well.  And  just  as  in  the  linear  theory,  these  incremental 
relations  are  exact  for  continuous,  piecewise-linear  strain  histories  (with  a  possible  jump 
discontinuity  at  t  =  0  only),  provided  that  any  jump  discontinuities  in  the  strain  rate  E  occur  at 
one  of  the  discrete  times 

We  may  use  the  conclusions  in  the  preceding  paragraph  to  derive  exact  relations  for  S(t)  for 
simple  piecewise-linear  strain  histories.  The  results  are  similar  to  those  in  section  3.4.5,  except 
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that  any  terms  there  involving  //^  require  appropriate  modification  to  account  for  the  fact  that  the 
equilibrium  elastic  response  function  Goo  is  no  longer  given  by  G 00(E)  =  2/ioo E.  For  example, 
for  the  constant  strain-rate  case  with  a  possible  initial  jump  in  strain  (cf.  equations  101-102),  the 
relation  104  for  S(t)  needs  to  be  changed  to 


S(t)  —  2^oo(E(t))  + 


'  N 

£ 

n=l 


2fMnTn(l-e  t,Tn ) 


E 1  + 


N 


0-t/Tn 


n=  1 


E{  0) 


(202) 


4.4.2  One-Parameter  Family  of  Functions:  Part  I 

Next,  we  consider  a  stress  relaxation  function  G  that  is  intermediate  in  generality  between  the  one 
in  the  the  previous  section,  equation  193,  and  the  stress  relaxation  function  for  the  general  Prony 
series  158  with  the  functions  unrelated.  While  the  flexibility  inherent  in  the  latter  case  allows 
for  a  wider  range  of  materials  to  be  approximated  than  would  be  the  case,  say,  with  the  separable 
stress  relaxation  function  in  equation  157  or  the  slightly  more  general  stress  relaxation  function  in 
equation  193,  when  it  comes  to  calibrating  the  model  for  a  particular  material  we  are  eventually 
forced  to  choose  specific  functional  forms  for  the  „  or  \Pn  (recall  that  \F.„(C)  :=  fly, (IE)).  This 
task  is  simplified  if  each  of  these  functions  is  assumed  to  belong  to  some  specific  subclass  of  the 
set  of  all  isotropic  functions  satisfying  equation  160,  and  in  particular  if  the  functions  in  this 
subclass  are  generated  by  a  single  real  parameter. 

Thus  let  and  $  be  isotropic  functions  satisfying 

*(u,C)  =  &(u,E)  =  E  +  0^(E2),  (203) 


where  u  is  a  non-dimensional  real  parameter.  And  for  each  for  n  —  1, . . . ,  N,  assume  that 
<E *n(E)  is  given  by  <&(u i,E)  and  \Pn(C)  by  \I> (u>,  C)  for  some  particular  choice  of  ui  which  may 
depend  on  n,  say  u  =  un.  Then  for  n  =  1, . . . ,  we  have 


®n(E)  =  &(u 0n,E) 

*n(C)  =  V((Jn,C)  , 


(204) 


and  by  equation  203  the  condition  160  is  satisfied.  We  do  not  necessarily  impose  the 
assumption  204  for  n  =  00,  although  as  a  special  case  we  might  do  so,  as  in  the  example  that 
follows.  Note  that  if  condition  204  is  not  imposed  for  n  =  00,  then  it  generally  will  not  hold  for 
n  =  0  either,  since  <F0  and  4? 00  arc  related  by  equation  167. 
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As  an  illustration  of  the  above,  we  take 


c)  =  |  (/  -  C -1)  +  (CT1  -  C~2) ,  0  <  u;  <  1 .  (205) 

Motivation  for  this  choice  is  discussed  below.  By  equation  4,  we  have 

C~2  =  (CT1)2  =  /  -  4E  +  0(E2) ,  (206) 

which  together  with  equation  4  implies  that 

1(1- C -1)  &  i(C_1-C-2)  =  £  +  0(£2).  (207) 

Then  on  substituting  equation  207  into  equation  205,  we  see  that  this  choice  for  \F(u;,  C)  indeed 
satisfies  equation  203.  We  will  assume  that  equation  204  holds  not  only  for  n  —  1 , ,N  but 
also  for  n  =  oo.  And  each  of  the  un  in  equation  204  must  satisfy  the  inequality  in  equation  205: 

0  <  oon  <  1  (208) 


forn  =  1, . . . ,  N,  oo.  Then  on  using  equations  204-206,  208,  167,  and  166,  we  find  that  the 
normalized  instantaneous  elastic  response  function  is  also  of  the  form  205.  More  precisely, 


3>0(E)  =  *0(C)  =  *(u0,C) 


(209) 


for 


o;0 


/A 

A*o 


N 

E 

n= 1 


/A 

Ao 


-Ur 


(210) 


and  LJ()  satisfies  the  inequality  0  <  uq  <  1.  Thus  equations  208  and  204  hold  for  n  —  0  also. 


From  the  relations  above,  equation  165,  and  equation  163,  we  see  that  the  instantaneous  and 
equilibrium  elastic  response  functions  are  given  by 


G0(E)  =  G0{C)  =  2  fio^(uo,C)  (211) 

and 

Goc(E)  =  Goo(C )  =  2/x00\Er(u ;oo,  C ) .  (212) 


This  example  is  motivated  by  the  Mooney-Rivlin  model  for  incompressible  elastic  materials.  We 
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give  a  brief  discussion  of  that  model  in  section  4.4.3  and  then  return  to  the  discussion  of  the 
example  above. 

4.4.3  The  Mooney-Rivlin  Elastic  Model 

The  Mooney-Rivlin  constitutive  model  for  an  isotropic,  incompressible,  elastic  solid  is  given  by14 

<r  =  pi  +  p+B  —  p~  B^1 ,  (213) 

where  p  is  indeterminate  and  the  coefficients  /i+  and  p~  are  elastic  moduli  satisfying 

/i+>0,  p~>0.  (214) 

B  is  the  left  Cauchy-Green  deformation  tensor: 

B  :=  FFt  =  V2  =  RCRt  ;  (215) 

and 

B~x  —  F~TF~l  —  V~2  .  (216) 


On  taking  the  small  strain  limit  of  equation  213,  one  finds  that  the  moduli  p+  and  p  are  related 
to  the  linear  elastic  shear  modulus  p  by 


/i  fl  -  /i  • 


(217) 


If  we  define  the  non-dimensional  parameter  u  by 


u  = 


IF 


P+  T  /r  P 


_  F_  _  x 


(218) 


then 


so 


p 

uj  >  0  and  0  <  —  =  1  —  u , 

p 


p+  =  pco,  p  —p{l  —  uj)  and  0  <  to  <  1 . 


(219) 


(220) 


14Cf.  Truesdell  and  Noll  (2,  §95).  The  indeterminate  scalar  p  is  sometimes  referred  to  as  the  “pressure”  and  often 
denoted  by  p,  although  it  is  not  equal  to  the  pressure  p  as  we  have  defined  it  in  equation  10i.  The  latter  is  given  by 

p  =  —p  —  |(/r+tr  B  —  p~  tr  B^1). 
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Thus  we  may  write  the  Mooney-Rivlin  constitutive  relation  213  as 


cr  =  —pi  +  p  [ujB  —  (1  —  c o)B  ,  0  <  u  <  1 .  (221) 

The  special  case  u  =  1  (equivalently,  p~  =  0)  is  called  a  neo-Hookean  material. 

Next,  note  that  by  equations  215,  216,  and  4i,  we  have 

F-1F-t  =  C~ \  F  lBF  T  =  I ,  F  lB  lF  T  =  C~2  .  (222) 


Then  by  equations  213,  221,  222,  and  7±  and  the  fact  that  J  =  1  for  an  incompressible  material, 
we  see  that  the  Mooney-Rivlin  constitutive  relation  is  equivalent  to 


s  =  -pc-1  +  p+i  -  p-c-2 
=  -p  cr 1  +  p  [ujI  -  (1  -  c o)C-2]  . 

This  may  be  re-written  in  the  equivalent  forms 


(223) 


5  =  -p* c -l  +  fj+  (/  -  cr x)  +  p-  (cr 1  -  cr2) 

=  -p^C^1  +  2yU\I >(cu,  C) , 


(224) 


where  p„  is  also  indeterminate  but  different  than  p,  and  ^(cu,  C )  is  given  by  equation  205. 


Though  the  relations  224  are  slightly  more  complicated  than  the  equivalent  relations  223,  they  are 
more  useful  for  our  purposes.  Indeed,  as  noted  in  the  previous  section,  the  function  T'  (jj.  C )  in 
equation  2242  satisfies  the  requirement  203.  On  the  other  hand,  if  we  rewrite  equation  2232  as 


S  =  -pc-1  +  2p-  2[cul- (l-cu)cr2]  , 


(225) 


we  find  that 

\[ul-  (1  ~u)C~2]  ±  E  +  0{E2).  (226) 


4.4.4  One-Parameter  Family  of  Functions:  Part  II 

Now  return  to  the  discussion  in  section  4.4.2.  When  the  assumptions  204  hold  for 
n  —  1, . . . ,  IV,  oo,  the  Prony  series  approximation  158  to  the  stress  relaxation  function  G (E.  t ), 
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when  expressed  in  terms  of  C,  takes  the  form 


N 

G(C,t)  =  2/u00^(o ;oc,C)  +  e-t/Tn2nn®{ujn,C) .  (227) 

n=  1 

As  mentioned  previously,  the  Mooney-Rivlin  model  provides  the  motivation  for  choosing  the 
one-parameter  family  of  functions  \P  (oj,  C)  to  have  the  particular  form  in  equation  205  .  Each 
term  in  this  Prony  series  involves  (the  determinate  part  of)  a  particular  Mooney-Rivlin 
constitutive  function,  as  is  clear  from  equation  2242. 

In  particular,  on  comparing  equations  211-212  with  equation  2242,  we  see  that  the  instantaneous 
and  equilibrium  elastic  response  functions  for  the  viscoelastic  part  S  of  the  model  are  given  by 
(different)  Mooney-Rivlin  constitutive  functions.  If  we  recall  the  decomposition  17, 

S  =  -JpC~l  +  S ,  (228) 

and  the  constitutive  assumptions  for  p  in  section  3.2.2,  and  compare  with  equation  2242,  we  see 
that  the  instantaneous  and  equilibrium  elastic  response  functions  for  the  total  2nd  Piola-Kirchhoff 
stress  tensor  S  are  given  by  (different)  compressible  versions  of  the  Mooney-Rivlin  model. 

Next,  we  examine  more  closely  the  form  of  the  stress  relaxation  function  in  equation  227  for 
arbitrary  t  when  \P (c u,  C )  is  given  by  equation  205.  For  this  choice,  the  term  2 //,, <P (ujni  C)  in 
equation  227  can  be  written  as 

2/^(u;„,  C)  =  rt(I  -  cr1)  +  Vn(C~x  -  C-2) ,  (229) 

where,  analogous  to  equation  220,  we  have  (for  n  —  1, . . . ,  N,  oo) 

P'ra  ■  ^ n  0  ,  [ln  .  fln  (1  U) n )  A  0  .  (230) 

Then  on  substituting  equation  229  into  equation  227,  we  see  that  the  stress  relaxation  function  is 
given  by 

G(C,  t )  =  -  C -1)  +  a^XC-1  -  C-2) ,  (231) 

where  the  functions  /i+  and  jr  are  given  by  the  Prony  series 

N 

■■=  nt,  +  YVnZ~t,Tn  ■  (232) 

n=  1 
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On  comparing  equation  231  with  equation  2241;  we  see  that  at  each  instant  t  the  function 
C  i — >  G(C ,  t)  is  (the  determinate  part  of)  a  particular  Mooney-Rivlin  constitutive  function.  It 
follows  that  the  total  stress  in  a  stress  relaxation  test  is  given  by  a  compressible  and 
“time-dependent”  version  of  the  Mooney-Rivlin  model. 

On  setting  t  =  oo  and  t  —  0  in  equation  232,  we  obtain 

V±(00)=  vt*  (233) 


^(0)  =  /4,  +  Vn  =■  Vo  •  (234) 

n=  1 

Then  on  using  equation  234,  the  definition  166  of  /r0,  and  the  definition  210  of  u0,  we  find  that 
equation  230  also  holds  for  n  —  0.  And  on  setting  t  —  0  and  t  =  oo  in  equation  231  and  using 
equations  234,  233,  230  (for  n  —  0  and  n  =  oo),  and  equation  205,  we  recover  equations  211  and 
212. 


The  stress  relaxation  function  in  equations  230-232  represents  a  substantial  simplification  of  the 
general  stress  relaxation  function  G(E,  s )  =  G(C .  s )  in  the  Pipkin-Rogers  model  as  well  as  a 
substantial  simplification  of  the  general  Prony  series  approximation  158  to  it.  However,  use  of 
equations  230-232  does  not  simplify  the  form  of  the  approximate  incremental  relations  in 
section  4.3.4.  On  the  other  hand,  the  evaluation  of  £n(tk)  at  each  of  the  discrete  times  A  is  fairly 
simple  once  C~l(tk)  has  been  determined.  Indeed,  by  equations  185,  204,  and  205, 

£n(tk)  =  *(ujn,C(tk)) 

/i  \  i  _  (235) 

=  T1  +  G  “  w")  C"(4)  -  • 


As  observed  in  section  4.3.1,  the  moduli  ji^  and  //„  and  the  relaxation  times  rn  for  the  general 
Prony  series  approximation  158  (and  hence,  in  particular,  for  the  special  approximation 
considered  here)  can  be  determined  from  the  viscoelastic  response  to  small  strain  deformations, 
for  which  the  linear  theory  in  section  3.  is  adequate.  For  the  particular  model  considered  in  the 
present  section,  the  departure  from  the  linear  theory  is  completely  characterized  by  the  additional 
material  parameters  aq, . . . ,  u>n,  tUoo.  Thus,  one  could  determine  the  moduli  //,,  and  relaxation 
times  rn  from  small  strain  tests  first,  and  then  optimize  the  choice  of  the  c on  to  fit  the  large  strain 
response. 
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